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ABSTRACT 

The  principal  focus  of  Automatic  Object  Recognition  (AOR)  involves  the  generation  of  appropriate  algorithms 
to  process  the  output  of  multi-spectral  sensor  arrays.  Given  the  high  dimensionality  that  characterizes  the 
signatures  of  targets  of  interest,  it  is  normally  impossible  to  satisfy  the  need  for  raw  signature  data  by  means  of 
measurement  records  a’.one.  Individual  sensor  characteristics  in  conjunction  with  aspect-angle  dependence,  target 
and  background  conjuration  (singly  and  in  synergism)  and  multi-spectral  tradeoffs  inexorably  lead  to  a 
requirement  for  predictive  signature  modeling  methods.  By  means  of  this  stratagem,  a  measured  signature  data 
base  can  be  leveraged  significantly,  improving  the  fidelity  of  the  overall  simulation. 

Irrespective  of  the  specific  representation  used  for  a  three-dimensional  geometry  and  material  database,  rarely 
does  a  predictive  signature  application  code  read  that  database  directly.  Rather,  a  specific  interrogation  method  is 
used  to  pass  particular  geometric  and  material  attributes  to  the  application  code.  Clearly  the  nature  of  the  physics 
employed  in  the  application  is  both  enabled  and  constrained  by  the  form  of  the  interrogation  process  used. 

In  this  paper,  several  examples  of  predictive  radar  codes  are  given,  illustrating  several  strikingly  different  ways 
of  linking  geometry  to  applications.  Following  those  examples  the  interface  methods  known  to  the  authors  will  be 
described  While  many  of  the  techniques  have  already  been  implemented,  some  are  currently  in  development.  In 
addition,  the  utility  of  various  techniques  will  be  related  to  particular  application  codes. 

1.  INTRODUCTION 

The  principal  focus  of  Automatic  Object  Recognition  (AOR)  involves  the  generation  of  appropriate  algorithms 
to  process  the  output  of  what  are  often  multi-spectral  sensor  arrays.  The  generation  and  testing  of  such 
algorithms  necessarily  require  the  full  panoply  of  desired  object  (target)  signatures  as  well  as  undesired  and 
unavoidable  (clutter)  signatures.  Even  ignoring  the  sensor  response  characteristics  (time  and  space  resolution, 
polarization,  noise,  etc.),  both  the  target  and  clutter  signature  sets  are  unbounded  even  for  a  single  frequency,  let 
alone  in  multi-spectral  domains.  The  unboundedness  is  strictly  true  for  these  signature  sets  and  most  probably  so 
even  for  derived  statistical  measures. 

Given,  therefore,  the  requirements  of  a  typical  AOR  project  for  copious  signature  files  for  algorithm  generation 
and  testing,  it  is  normally  impossible  to  satisfy  the  need  by  means  of  measurement  records  alone.  The  single- 
frequency  sensor  characteristics  noted  above  in  conjunction  with  aspect  angle  dependence,  target  and  background 
configuration  (singly  and  in  synergism)  and  multi-spectral  tradeoffs,  inexorably  lead  to  a  requirement  for  predictive 


+  Presented  it'  the  Advanced'  Institute  rrogrum  on  Automatic  OBjccfc  Recognition,  sponsored  by  the  bocietj’  of  Photo-Optica! 
Instrumentation  Engineers  (SPIE),  Cocoa  Beach,  FL,  April  21-23,  1990. 
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modeling  methods.  By  means  of  this  stratagem,  a  measured  signature  data  base  can  be  leveraged  significantly, 
improving  the  fidelity  of  the  overall  simulation. 

For  more  than  twenty  years,  the  Ballistic  Research  Laboratory  (BRL)  has  been  utilizing  solid  geometric 
modeling  methods  to  support  vulnerability,  lethality  and  neutron  transport  studies  of  military  targets.  In  such 
item-level  studies,  target  geometry  and  material  information  are  passed  to  various  application  codes  to  derive 
certain  measures-of-performance.  Building  on  the  general  paradigm,  workers  at  BRL  and  elsewhere  have  extended 
the  general  techniques  to  support  many  predictive  signature  models4'6  including  optical,  millimetre  wave 
(MMW)/  Infra-Red  (IR),  magnetic  and  X-ray  models. 

It  is  important  to  note  that  this  type  of  analysis  must  generally  be  supported  by  a  solid  geometric  model.  A 
solid  model7  is  a  computer  description  of  closed,  solid,  three-dimensional  shapes  represented  by  an  analytical 
framework  within  which  the  three-dimensional  material  can  be  completely  and  unambiguously  defined.  Two  major 
families  of  solid  model  representations  exist,  each  with  several  unique  advantages.  The  first  is  the  Combinatorial 
Solid  Geometry  Representation  (CSG-Rep).8  Solid  models  of  this  type  are  expressed  as  boolean  combinations  of 
primitive  solids  which  are  geometric  entities  described  by  some  set  of  parameters  and  occupying  a  fixed  volume  in 
space.  The  second  is  the  Boundary  Representation  (B-Rep),  of  which  there  are  two  sub-types:  (l)  explicit,  where 
each  solid  is  described  by  an  explicit  enumeration  of  the  extent  of  the  surface  of  the  solid;  and  (2)  implicit,  where 
the  surface  of  the  solid  is  described  by  an  analytic  function  such  as  a  Coons  patch,  Bezier  patch,  B-spline,  etc. 
Hybrid  systems  such  as  the  BRL-CAD  Package9  also  exist. 

Irrespective  of  the  specific  representation  used  for  a  three-dimensional  geometry  and  material  data  base,  rarely 
does  an  application  code  read  that  database  directly.  Rather,  a  specific  interrogation  method  is  invoked  to  pass 
particular  geometric  and  material  attributes  from  a  source  or  reference  file  to  the  application  code.  Clearly  the 
nature  of  the  physics  employed  in  the  application  is  both  enabled  and  constrained  by  the  form  of  the  interrogation 
process  used. 

In  this  paper,  several  examples  of  predictive  radar  codes  are  given,  illustrating  several  strikingly  different  ways 
of  linking  geometry  to  applications.  Following  those  examples  the  interface  methods  known  to  the  authors  will  be 
described.  Many  of  the  techniques  have  already  been  implemented;  some  are  currently  in  development.  In 
addition,  the  utility  of  various  techniques  will  be  related  to  particular  application  codes. 


2.  CASE  STUDIES  IN  PREDICTIVE  RADAR  MODELING 
2.1  Examples  of  Synthetic  Aperture  Radar  Imagery 

Early  radars  were  designed  to  estimate  target  range  and  closing  rates.  For  these  systems,  the  main  parameter 
of  interest  for  a  given  target  was  Radar  Cross  Section  (RCS).  The  RCS  figure  represents  the  efficiency  with  which 
radar  waves  are  scattered  back  to  the  receiver.  Certain  modern  radars,  when  placed  on  moving  platforms  such  as 
aircraft,  can  be  used  to  form  two-dimensional  images  of  targets.  Radar  imagery  of  this  class  is  called  Synthetic 
Aperture  Radar  (SAR).  A  description  of  an  Armored  Fighting  Vehicle  (AFV)  has  been  analyzed  with  a  SAR 
program10  (to  be  described  below)  and  the  results  are  shown  in  Figs.  1  and  2. 

Figure  1,  left-hand  side,  shows  the  AFV  as  seen  by  the  SAR  radar  from  a  (35°,  30°)  (Azimuth,  Elevation) 
orientation.  A  horizontal  flight  path  (left  to  right)  is  assumed.  The  properties  of  SAR  processing  are  such  that 
following  signal  detection  and  manipulation  an  image  is  derived  which  resolves  the  target  in  range  and  cross-range 
(along  the  flight  path)  but  not  in  the  remaining  orthogonal  direction.  Thus  the  final  SAR  image  orientation  is 
similar  to  the  optical  rendering  shown  in  the  right-h  ,nd  side  of  Fig.  1. 


t  Recently,  an  initiative  was  made  to  consolidate  Radar  Cross  Section  (RCS)  code  development  sponsored  by  the  Tri-Services  and  NASA. 
The  Electromagnetic  Code  Consortium  has  chosen  the  geometry  tools  developed  by  the  BRL  as  the  basis  for  radar  codes  to  be  sponsored  by 
this  group.  See  Ref.  6. 
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Figure  x.  On  the  left  an  optical  image  of  an  AFV  i'  astrating  the  radar  view  of  the  target  (35°,  30°) 
for  a  synthetic  aperture  radar  (SAR)  simulation.  The  SAR  is  modeled  as  moving  in  the  azimuthal 
direction  (elevation  and  range  constant}.  On  the  right,  complementary  optical  image  to  that  shown 
on  the  left,  the  apparent  aspect  is  (Xl5°,  80°^  and  is  i.  aggestive  of  the  SAR  reconstructive  process 
when  range  is  plotted  against  crosa-range  as  in  Fig.  2.  ^Frorn  Ref.  10) 


Figure  2,  Computed  SAR  image  for  target  resolution  of  10  inches.  Left-hand  image 
Vertical/  Vertical  (co-polarization  transmit/  receive)  modes;  right-hand  image  is 
Horizontal  (cross-polarizati  m  transmit/  receive)  modes.  (From  Ref.  10) 


is  for  10  Ghz, 
for  Vertical/ 
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A  pair  of  SAR  images  for  the  AFV  is  shown  in  Fig.  2.  Both  images  are  computed  for  a  target  resolution  of  10 
inches.  The  left,  image  is  the  result  of  a  transmit  Vertical,  receive  Vertical  polarization  mode;  the  right,  for  a 
Vertical/Horizontal  polarization  mode.  In  each  of  these  images,  the  radar  signal  is  propagating  from  top  down. 
Range  information  is  plotted  along  the  ordinate  and  cross-range  data  along  the  abscissa. 

2.2  The  SRIM  Code  of  ERIM 


Simulated  Radar  IMagery  (SRIM)  is  a  high-frequency  predictive  radar  code  developed  at  the  Environmental 
Research  Institute  of  Michigan  (ERIM).11  Early  versions  of  the  code  used  a  geometry  representation  scheme  for 
which  no  convenient  editor  existed.  In  1983,  ERIM  linked  its  electromagnetics  section  of  SRIM-to  BRL’s  CSG-Rep 
geometric  modeling  capability.  BRL  geometry  is  described  in  terms  of  Boolean  combinations  (intersections, 
unions,  differences)  of  primitive  objects  such  as  ellipsoids,  c;  :ders,  and  tori.  Ray-tracing  is  used  to  extract 
geometric  information  from  the  database.  The  ERIM-BRL  link  provides  an  easy  means  for  generating  geometry 
descriptions  through  an  excellent  graphics  editor  and  a  natural  geometry  interrogation  technique  for  determining 
ray  path  information.  SRIM  follows  the  Geometric  Optics  (GO)  paths  of  rays  as  they  hit  and  reflect  from  the 
target  surfaces.  For  each  ray  emanating  from  the  radar  a  piecewise-linear  path  is  traced  through  the  geometry  ?nd 
a  history  is  recorded  of  hit-point  coordinates,  normal  vectors,  principal  curvatures,  principal  curvature  directions 
and  a  flag  designating  whether  or  not  line-of-sight  exists  to  the  radar.  This  geometric  information  is  thciupasscd  to 
the  electromagnetics  portion  of  the  model.  At  each  hit  point  along  the  ray  path,  the  incoming  illumination 
wavefront  is  resolved  into  parallel  and  perpendicular  components,  a  complex  Physical  Optics  (PO)  field 
contribution  at  the  radar  is  calculated,  taking  into  account  the  wavelength  of  the  system  and  this  contribution  is 
then  resolved  into  its  vertical  and  horizontal  polarization  components.  SRIM  is  thus  able  to  determine  the  PO 
field  contributions  for  not  only  first-surface  illumination,  but  also  for  contributions  due  to  multiple  bounces  along 
the  specular  directions.  If  a  total  radar  cross  section  is  desired,  this  complex  return  is  summed  into  a  total  field 
variable.  If  a  SAR  image  is  desired,  the  location  at  which  this  contribution  would  appear  in  an  image  is  calculated 
and  added  to  the  appropriate  range/cross-range  bin.  Image  parameters  such  as  resolution,  beam  weighting  and 
pixel  size  are  specified  in  a  separate  file. 

If  a  vehicle  on  a  ground  plane  is  the  object  of  interest,  SRIM  can  generate  a  clutter  model  for  the  ground 
backscatter  and  will  properly  represent  the  shadowing  of  the  ground  plane  by  the  vehicle.  In  fact,  all  shadowing 
effects  are  correctly  represented  as  a  natural  result  of  the  ray-tracing  paradigm  used  in  this  simulation. 

2.3  The  SARSIM  Code  of  Northrop 


SARSIM  is  another  high-frequency  synthetic  aperture  radar  simulation  code.  It  was  developed  at  the  Northrop 
Research  and  Technology  Center."  In  this  model,  the  radar  images  of  the  target  and  surrounding  background  are 
computed  separately.  In-a  final  step  the  two  images  are  combined  via  processing  which  reflects  appropriate  noise, 
speckle  and  system-response  values.  For  the  geometric  calculations  an  underlying  target  representation  of 
triangular  plates  is  employed.  For  a  given  SAR  resolution,  each  of  the  triangular  plates  in  the  target  description  is 
subdivided  into  panels  whose  contributions  fall  into  a  single  SAR  resolution  bin.  For  each  of  these  panels  in  the 
target  description,  a  PO  calculation  is  carried  out  to  determine  the  contribution  to  the  corresponding  SAR 
resolution  bin.  This  is  also  done  for  a  Physical  Theory  of  Diffraction  calculation  on  the  panel  edges.  As  of  1987 
only  the  specular  (first  surface)  and  diffraction  contributions  were  considered  in  this  model.  A  simplified-shadowing 
algorithm  eliminates  any  panels  partially  obscured  by  other  panels.  Multiple  reflections  are  not  considered  yet. 

2.4  The  TRACK  Code  of  GTRI 


Georgia  Tech  Research  Institute  (GTRI)  has  developed  a  radar  prediction  tool  called  TRACK.  The 
geometry  used  by  TRACK  is  based  on  a  subset  of  objects  supported  by  the  GTRI  MAX  geometry  editor.  A  hybrid 
representation  of  simple  shapes  such  as  ellipsoids,  plates,  frusta  and  special  radar-specific  constructs  such  as 
dihedrals  and  trihedrals  are  accepted  by  the  TRACK  code.  TRACK  has  a  series  of  subroutines  referred  to  as 
CROSS'which  calculate  the  field  contribution  of'the  individual  scatteieis  in  a  MAX  geometry  file.  CROSS  uses  PO 


to  calculate  the  fields  for  polygonal  flat  plates,  right  circular  frusta,  general  ellipsoids,  ogives  and  toroids.  A 
combination  of  Geometric  and  Physical  Optics  is  use  to  calculate  the  fields  for  dihedral  and  trihedral  shapes. 
CROSS  also  predicts  diffraction  from  straight  and  curved  edges  by  methods  of  equivalent  currents  to  avoid 
problems  associated  with  caustics.  The  returns  from  the  individual  scatterers  can  be  summed  coherently, 
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noncoherently,  or  retained  in  complex  form  for  postprocessing.  Shadowing  is  handled  by  ignoring  all  plates  that 
face  away  from  the  radar  and  by  excluding  from  the  electromagnetics  calculations  any  scatterer  that  is  completely 
obscured  by  other  geometric  elements.  When  an  element  is  partially  covered,  the  full  energy  of  that  element  is 
returned,  which  can  sometimes  be  a  source  of  error.  For  post-processing,, a  file  with  scatterer  types,  location,  size 
and  orientation  is  generated.  This  can  then  be  used  in  numerous  scenarios,  such  as  SAR,  doppler  radar,  and  ISAR 
(Inverse  SAR). 

2.5  Discussion  of  Approaches 

Each  of  the  aforementioned  high-frequency  scattering  codes  required  compromise.  The  choice  of  geometric 
representation  and  interrogation  method  led  to  advantages  and  shortcomings.  Choosing  a  particular  geometric 
approach  often  requires  approximations  in  the  implementation  of  the  physics  of  an  analysis.  Likewise,  a  particular 
formulation  of  a  physical  analysis  can  force  the  acceptance  of  a  simplified  and  less  general  geometry  representation 
for  the  sake  of  getting  some  portion  of  the  problem  solved. 

SRIM  uses  the  information  gathered  by  ray-tracing  to  determine  a  scattered  field.  Fields  are  calculated  locally 
so  that  the  ray-trace  sampling  density  can  be  chosen  small  enough  to  guarantee  that  each  contribution  to  the 
scattered  field  will  lie  in  a  unique  resolution  cell,  a  feature  not  guaranteed  by  the  GTRI  model.  While  the  ray¬ 
tracing  technique  is  computationally  intensive,  the  automatic  shadowing  and  multiple  bounce  calculations  are 
attractive  features.  If  two  modeled  surfaces  are  at  a  right  angle  to  each  other,  the  SRIM  code  will  calculate  a 
scattered  field  that  is  appropriate  for  a  dihedral  without  requiring  the  explicit  designation  of  a  dihedral  element  in 
the  geometry  file.  This  approach  does  have  the  severe  cost  of  requiring  ray-tracing  densities  which  are  frequency 
dependent,  so  that  the  number  of  rays  which  must  be  cast  at  the  target  description  for  a  94  GHz  radar  calculation 
is  on  the  order  of  100  times  the  number  of  rays  needed  at  10  GHz.  Originally  the  code  simply  point-sampled  the 
field  contribution  at  a  location,  but  this  approach  suffered  from  aliasing  problems  due  to  under  sampling  high 
spatial  frequency  (fine  detail)  geometric  elements.  Improvements  have  been  made  so  that  the  field  contribution  is 
integrated  over  a  planar  approximation  at  the  scattering  surface,  so  that  field  calculations  are  less  sensitive  to  ray- 
trace  sampling  densities.  Even  so,  with  a  ray-tracing  approach,  one  can  still  run  the  risk  of  not  adequately 
sampling  the  fine  detail  in  the  geometric  model,  although  increasing  sampling  densities  with  ray-tracing  makes  the 
inter-ray  sampling  distances  arbitrarily  small.  It  is  possible  to  place  an  upper  bound  on  the  size  of  geometric 
objects  that  will  be  under  sampled;  with  knowledge  about  the  distribution  of  feature  sizes  in  the  geometry,  the  ray 
spacing  can  be  chosen  so  as  to  make  it  statistically  unlikely  to  miss  any  details. 

Only  PO  fields  are  calculated  for  SRIM.  Diffraction  effects  occurring  at  edges  are  not  calculated  because  the 
ray-tracing  interrogation  method  does  not  provide  the  required  edge  information.  Some  of  the  representational 
schemes  discussed  later  will  show  how  this  might  be  accomplished  with  an  alternate  form  for  extracting 
information  from  the  model  geometry. 

One  advantage  of  the  SARSIM  approach  to  field  calculations  is  that  integration  takes  place  over  each  facet  or 
panel  in  a  target  description  and  all  of  the  modeled  geometry  is  guaranteed  to  be  included  in  the  calculation.  SRIM 
can  not  guarantee  this.  Likewise,  if  one  has  adopted  a  particular  resolution' cell  size  and  subdivided  the  triangular 
plates  into  panels  appropriate  for  the  resolution  cell  size,  there  will  be  no  increase  in  the  number  of  panels  for  an 
increase  in  the  frequency  of  the  radar.  Gn  the  other  hand,  in  modeling  the  target  with  flat  plates,  one  must 
approximate  any  curved  surfaces  and  thereby  introduce  some  uncertainty  about  the  fidelity  of  the  geometry  itself, 
since  flat  plate  representations  of  curved  surfaces  introduce  edges  where  no  real  edges  exist.  Thus,  each  edge  must 
be  marked  by  the  geometry  modeler,  to  indicate  if  the  edge  is  an  actual  edge,  or  a  modeling  artifact.  For  radar 
diffraction  calculations  this  criterion  must  be  used  to  trigger  the  inclusion  or  exclusion  of  edges. 

Multiple  reflections  are  not  modeled  in  SARSIM,  so  any  radar  signature  dominated  by  multiple  bounces  of 
energy  within  the  target  will  be  poorly  modeled.  Many  ground  vehicles,  especially  tracked  vehicles,  exhibit 
dramatic  dihedral  and  trihedral  returns,  making  multiple-bounce  considerations  a  primary  concern.  The  use  of 
piecewise-linear  target  representations  does  lend  itseif  to  polygonal  projection  for  determining  multiple  path 
interaction,  thereby  avoiding  ray-tracing.  However,  such  an  approach  introduces  its  own  set  of  difficulties.  These 
techniques  require  significant  computation  to  find  all  the  facet-to-facet  “form  factors”.  For  a  geometry  with  N 
facets,  a  solution  to  the  global  illumination  problem  requires  a  system  of:  N  equations  in  N  unknowns  to  be  solved. 
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With  most  targets  of  interest  the  value  of  N  is  quite  large,  and  performing  Gaussian  elimination  on  the  N-by-N 
matrix  is  often  prohibitive. 

The  GTRI  radar  simulation  provides  highly  accurate  electromagnetic  calculations  for  most  of  its  primitive 
geometry  types.  Of  the  three  models  considered  here,  it  provides  the  most  detailed  calculations  for  certain  simple 
objects.  For  many  modeling  scenarios  this  approach  gives  excellent  results.  However,  the  technique  for  combining 
these  component  field  calculations  for  a  complex  geometric  object  fails  to  address  some  important  issues.  The  basic 
approach  of  modeling  complex  geometry  with  simple  objects  for  which  highly  detailed  electromagnetic  field 
calculation  are  known,  has  been  used  for  many  years.14  For  simple  shapes,  e.g.  rockets,  or  for  objects  where  a  few 
scatterers  with  wide  spatial  separation  dominate  the  returns,  this  has  been  a  highly  successful  technique.  When  the 
target  under  consideration  has  non-convex  complexity  such  as  seen  in  a  tank,  the  effects  of  multiple-reflection 
scattering  cannot  be  ignored.  Although  the  GTRI  model  has  dihedral  and  trihedral  geometric  primitives  to 
account  for  some  multiple  reflection  scattering,  the  modeler  must  provide  the  information  of  how  large  and  where 
these  objects  should  be.  Even  with  this  feature,  multiple  reflection  paths  may  escape  the  modeler’s  attention  or  be 
of  such  complexity  as  to  be  insufficiently  modeled  by  dihedrals  or  trihedrals.  To  deal  with  this  type  of  scattering 
with  a  general-purpose  code  that  does  not  require  operator  intervention,  it  would  be  necessary  to  appeal  to  another 
geometric  model  of  higher  resolution  to  find  the  multiple  reflection  paths. 

One  other  shortcoming  of  this  method  is  the  coarseness  of  the  criterion  for  determining  when  one  primitive 
obscures  or  shadows  another.  The  GTRI  method  does  not  sufficiently  subdivide  the  scattering  geometry,  so  that 
full-power  returns  to  the  radar  are  considered  to  exist  even  when  the  scattering  object  is  almost  completely 
obscured.  Only  when  a  primitive  is  totally  obscured  is  the  electromagnetic  scattering  omitted  from  the  field 
calculation.  Notice  that  even  if  it  were  geometrically  simple  to  determine  what  portion  of  a  given  solid  is 
obscured,  this  electromagnetics  calculation  method  would  not  provide  the  field  contribution  from  the  partially 
visible  portion.  Still  worse  problems  arise  for  imaging:  radar  simulations  if  the  scattering  object  is  geometrically 
larger  than  a  resolution  cell,  since  a  given  scatterer  has  only  one  field  value  and  that  field  value  cannot  be 
distributed  over  more  than  one  cell.  This  all-or-nothing  field  calculation  method  is  a  consequence  of  the  physics 
being  constrained  by  the  chosen  geometry  interface. 

A  geometric  representation  method  and  an  interface  of  that  method  to  the  scattering  calculations  have  been 
selected  in  each  of  the  three  codes  just  discussed.  These  choices  clearly  have  advantages  and  liabilities.  Any  code 
requiring  geometric  information  will  likewise  be  limited  in  some  areas  and  enhanced  in  others  by  the  choices  of  the 
geometric  representation  and  the  interface  of  that  geometry  to  the  analysis.  Each  analysis  code  described  so  far 
employs  only  one  of  these  geometry/analysis  code  links,  but  significant  benefits  might  be  gained  from  analysis 
codes  that  use  two  or  more  of  these  interfaces  simultaneously.  Examples  of  the  information  that  can  be  readily 
provided  from  various  geometry  representations  and  interfaces  will  be  discussed  in  the  following  sections. 

3.  GEOMETRIC  INTERFACES 

In  what  follows  a  number  of  methods  are  described  by  which  geometric  and  material  data  are  supplied  to 
applications  codes.  Where  possible,  the  interface  methods  developed  and  supported  by  the  BRL-CAD  geometry 
tools9,15  will  be  used  to  describe  and  illustrate  these  processes.  A  goal  of  BRL-CAD  has  been  to  provide  a  general 
and  open  analysis  environment  in  which  a  variety  of  interrogation  interfaces  are  supported  so  that  diverse 
application  codes  can  be  driven  from  a  single,  unified  geometric  model.16  The  APPENDIX  gives  a  brief  overview 
of  the  BRL-CAD  modeling  environment  and  some  strategies  which  have  guided  its  development. 

3.1  Ray-tracing 

Rays  begin  at  a  point  1^,  and  proceed  infinitely  in  a  given  direction  given  by  the  unit  normal  vector!^.  The 
direction  vector  or  direction  cosines  for  the  ray  (Dx,Dy,D,)  are  the  cosines  of  the  angle  between  the  ray  and  each  of 
the  Cartesian  axes.  Any  point "X  on  a  ray  may  be  expressed  as  a  linear  combination  ofd?  and  I?  by  the  formula 

t  =~P  +  k 

where  valid  values  fork  are  in  the  range  (  0,  oo]. 
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The  traditional  approach  to  ray-tracing  has  been  batch-oriented,  with  the  user  defining  a  set  of  “viewing 
angles”,  initiating  a  large  batch  job  to  compute  all  the  ray  intersections  and  then  post-processing  all  the  ray  data 
into  some  meaningful  form.  However,  the  major  drawback  of  this  approach  is  that  the  application  has  no 
immediate  control  over  ray  paths,  making  another  batch  run  necessary  for  each  level  of  reflection,  etc. 

In  order  to  be  successful,  applications  need:  (1)  interactive  control  of  ray  paths,  to  naturally  implement 
reflection,  refraction  and  fragmenting  into  multiple  subsidiary  rays  and  (2)  the  ability  to  fire  rays  in  arbitrary 
directions  from  arbitrary  points.  Nearly  all  non-batch  implementations  have  closely  coupled  a  specific  application 
(typically  a  model  of  illumination)  with  the  ray-tracing  code,  allowing  efficient  and  effective  control  of  the  ray 
patho.  The  most  flexible  approach  of  all  is  to  provide  the  ray-tracing  capability  through  a  general-purpose  library 
and  make  the  functionality  available  as  needed  to  any  application.  For  example,  the  decision  of  when  a  ray  should 
be  reflected,  transmitted,  or  absorbed  should  be  entirely  under  the  control  of  the  application  program.  A  set  of 
sample  ray  histories  that  mjght  result  from  an  application  exerting  such  interactive  control  is  shown  in  Fig.  3. 

3.1.1  RT  Library  Interface:  The  third  generation  ray-tracing  capability  in  the  BRL-CAD  Package  is  a  set  of 
library  routines  in  librt  to  allow  application  programs  to  intersect  rays  with  model  geometry.  There  are  two  parts 
to  the  interface:  “preparation”  routines  and  the  actual  ray-tracing  routine.  rt_dirbuild()  opens  the  database  file 
and  builds  the  in-core  database  table  of  contents.  rt_gettree()  adds  a  database  sub-tree  to  the  active  model 
space,  and  can  be  called  multiple  times  to  join  different  parts  of  the  database  together. 

To  compute  the  intersection  of  a  ray  with  the  geometry  in  the  active  model  space,  the  application  must  call 
rt_shootray()  once  for  each  ray.  Ray  behaviors  such  as  perspective,  reflection,  refraction,  etc,  are  entirely 
determined  by  the  applications  program  logic,  and  not  by  the  ray-tracing  library.  The  ray-path  specification 
determined  by  the  applications  program  is  passed  as  a  parameter  to  rt_shootray()  in  the  application  structure, 
which  contains  five  major  elements:  the  vector  a_ray.r_pt  (?)  which  is  the  starting  point  of  the  ray,  the  vector 
a_ray.r_dir  (?)  which  is  the  unit-length  direction  vector,  the  pointer  *a_hit()  to  an  application-provided  routine  to 
be  called  when  some  geometry  is  hit  by  the  ray,  the  pointer  *a_miss()  to  an  application-provided  routine  to  be 
called  when  the  ray  does  not  hit  any  geometry,  and  the  variable  a_onehit.  In  addition,  there  are  various  locations 
for  applications  to  store  state  information  such  as  recursion  level,  intermediate  color  values,  and  cumulative  ray 
distance. 

When  the  a_onehit  variable  is  set  to  zero,  the  ray  is  traced  through  the  entire  model.  Applications  such  as 
lighting  models  may  often  only  be  interested  in  the  first  object  hit;  in  this  case,  a_onehit  may  be  set  to  the  value 
one  to  stop  ray-tracing  as  soon  as  the  ray  has  intersected  at  least  one  piece  of  geometry.  Similarly,  if  only  the  first 
three  hits  are  required  (such  as  in  the  routine  that  refracts  light  through  glass),  then  a_onehit  may  be  given  the 
value  of  three.  Then,  at  most  three  hit  points  will  be  returned,  an  in-hit,  an  out-hit  and  a  subsequent  in-hit. 
When  only  a  limited  number  of  intersections  are  required,  the  use  of  this  flag  can  provide  a  significant  savings  in 
run-time. 

The  rt_shootray()  function  is  designed  for  full  recursion  so  that  the  application  provided  a_hit()/a_miss() 
routines  can  themselves  fire  additional  rays  by  recursively  calling  rt-shootr&yQ  before  deciding  their  own  return 
value.  In  addition,  the  function  rt_shootray()  is  fully  capable  of  operating  in  parallel  with  other  instances  of 
itself  in  the  same  address  space,  allowing  the  application  to  take  advantage  of  parallel  hardware  capabilities  where 
such  exist. 

3.1.2  Sample  RT  Application:  A  simple  application  program  that  fires  one  ray  at  a  model  and  prints  the 
result  is  included  below,  to  demonstrate  the  simplicity  of  the  interface  to  librt. 
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struct  application  ap; 
main()  { 

rt_dirbuild("model.g"); 

rt_gettree("car"); 

rt_prep(); 

VSET(  ap.a-point,  100,  0,  0  ); 

VSET(  ap.a_dir,  -1,0,0  ); 
ap.a_hit  =  &hit_geom; 
ap.a_miss  =  &miss_geom; 
rt_shootray(  <&ap  ); 

} 

hit_geom(app,  part) 
struct  application  *app; 
struct  partition  *part; 

{ 

printf("Hit  %s",  part->pt_forw->pt_regionp->reg_name); 

} 

miss_geom(){ 

printf("Missed"); 

} 

3.1.3  Ray  Intersection  Data:  If  a  given  ray  hits  some  model  geometry,  the  a_hit()  routine  is  called  and  is 
provided  a  pointer  to  the  head  of  a  doubly-linked  list  of  partition  structures.  Each  partition  structure  contains 
information  about  a  line  segment  along,  the  ray;  the  partition  has  both  an  “in”  (pt_inhit)  and  an  “out” 
(pt-outhit)  hit  point.  Each  hit  point  is  characterized  by  the  hit  distance  hit-dist,  which  is  the  distance  k  from 
the  starting  point  r_pt  along  the  ray  to  the  hit  point.  The  linked  list  of  partition  structures  is  sorted  by 
ascending  values  of  hit-dist.  As  a  result  of  this  definition,  the  “line-of-sight”  distance  between  any  two  hit  points 
can  be  determined  simply  by  subtracting  the  two  hit-dist  values.  This  will  give  the  distance  between  the  hit 
points,  in  millimeters. 

If  the  variable  a_onehit  was  set  non-zero,  then  only  the  first  a_onehit  hit  points  along  the  partition  list  are 
guaranteed  to  be  correct;  any  additional  hit  points  provided  should  be  ignored.  This  is  usually  important  only 
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when  a_onehit  was  set  to  an  odd  number;  in  this  case  the  value  of  pt-outhit  in  the  last  partition  structure  may 
not  be  accurate  and  should  be  ignored. 

If  the  actual  3-space  coordinates  of  the  hit  point  are  required,  they  can  be  computed  into  the  hit-point 
element  with  the  C-language  version  of  =  ^+  k  *1? : 

VJOINl(  hitp-> hit-point,  rayp->r_pt,  hitp->hit_dist,  rayp->r_dir ); 


3.1.4  Surface  Normals:  As  an  efficiency  measure  only  the  hit  distances  are  computed  when  a  ray  is  intersected 
with  the  model  geometry.  For  any  hit  point,  the  surface  normal  at  that  point  can  be  easily  acquired  by  executing 
the  C  macro: 

RT_HIT_NORM(  hitp,  stp,  rayp  ); 

In  addition  to  providing  the  unit-length  outward-pointing  surface  normal  in  struct  hit  element  hit-normal,  this 
macro  also  computes  the  3-space  coordinates  of  the  hit  point  in  struct  hit  element  hit-point. 

3.1.5  Gaussian  Curvature:  For  any  hit  point,  after  the  surface  normal  has  been  computed,  the  Gaussian 
surface  curvature  at  that  hit  point  can  be  acquired  by  executing  the  C  macro: 

RT_CURVE(  curvp,  hitp,  stp  ); 

At  the  hit  point,  there  exists  exactly  one  pair  of  orthogonal  directions  also  orthogonal  to  the  surface  normal  7$  for 
which  the  values  of  c  take  on  the  minimum  and  maximum  values  q  and  c2.  q  and  c2  are  the  inverse  radii  of 
curvature  and  |q  |<|c2|,  i.e.  q  is  the  most  nearly  flat  principle  curvature.  A  positive  curvature  indicates  that  the 
surface  bends  toward  the  (outward  pointing)  normal: vector^  at  the  hit  point.  A  curvature  structure  has  three 
elements,  the  unit  vector  crv_pd‘ir  (or  A)  pointing  in  the  direction  of  the  first  principle  curvature,  the  scalar 
crv_cl  (or  q)  gjving  the  curvature  in  the  first  principle  direction  and  the  scalar  crv_c2  (or  c2)  giving  the 
curvature  in  the  second  principle  direction  eT.  The  secondprinciple  direction  is  implied  and  can  be  found  by 
taking  the  cross  product  of  the  normal  with  crv_pdir,  i.e.,  B  =  X^. 

3.1.8  U-V  Mapping:  Each  primitive  solid  can  be  considered  to  be  bounded  by  one  or  more  regular  surfaces. 
Each  regular  surface  is  defined  as  the  locus  of  points  S(u,v)  depending  on  two  real  parameters  u  and  v  which  range 
from  0.0  to  1.0  inclusively.  These  parameters  form  the  coordinates  of  a  two-dimensional  Cartesian  u,v-plane.  A 
given  (u,v)  coordinate  will  appear  only  once  on  each  regular  surface,  but  in  objects  with  more  than  one  surface  that 
same  (u,v)  pair  may  appear  at  more  than  one  place  on  the  object.  The  (u,v)  coordinate  of  the  hit  point  is  returned 
in  uvcoord  structure  elements  uv_u  and  uv_v.  For  any  hit  point,  after  the  value  of  hit-point  has  been 
computed,  the  (u^v):coordinates  of  that  point  can  be  acquired  by  executing  the  C  macro: 

RT_HIT_UVCOORD(  ap,  stp,  hitp,  uvp  ); 

For  some  simple  optical  rendering  applications  ,  it  is  sometimes  desirable  to  create  a  mapping  between  the 
coordinate  system  on  the  surface  of  an  object  to  coordinates  on  a  plane.  This  is  generally  used  to  drive  simple,  two 
dimensional  texture  mapping  algorithms.  The  most  common  application  is  to  extract  a  “paint”  color  from  a 
rectangular  image  file  at  coordinates  (u,v),  and  apply  this  color  to  the  surface  of  an  object.  These  parameters  can 
also  be  used  to  simulate  the  effect  of  minor  surface  roughness  using  the  bump  mapping  technique.  Here,  the  u  and  V 
coordinates  index  into  a  rectangular  file  of  perturbation  angles;  the  surface  normal  returned  by  RT_HIT_NORM() 
is  then  modified  by  up  to  ±90  degrees  each  in  both  the  u  and  v  directions,  according  to  the  stored  perturbation. 

In  addition,  the  approximate  “beam  coverage”  of  the  ray  in  terms  of  the  parr  neters  (u,v)  is  returned  in  the 
structure  elements  uv_du  and  uv_dv.  These  approximate  values  are  based  upon  the  ray’s  initial  beam  radius 
(a_rbeam)  and  beam  divergence  per  millimeter  (a_diverge)  as  specified  in  the  application  structure.  These  delta-u 
and  delta-v  values  can  be  helpful  for  anti-aliasing  or  filtering  areas  of  the  original  texture  map  to  produce  an  “area 
sample”  value  for  the  hit  point. 
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3.2  Topological  Representation 

Some  predictive  radar  signature  codes,  such  as  the  TRACK  code  of  GTRI  discussed  above,  do  not  operate 
directly  on  a  solid  geometric  representation  of  an  object.  Instead,  they  rely  on  the  fact  that  large  radar  returns 
occur  primarily  due  to  the  existence  of  dihedral  and  trihedral  structures  in  the  object.  Rather  than  describing  a 
vehicle  simply  as  a  collection  of  these  topological  structures,  one  can  analyze  a  three  dimensional  solid  model  to 
locate  all  instances  of  the  topological  features  of  interest.  For  example,  the  software  could  locate  planar  face 
elements,  edges  where  two  locally  planar  elements  join  to  make  a  dihedral,  edges  where  three  locally  planar 
elements  join  to  make  a  trihedral,  etc.  Then  this  list  of  topological  features  becomes  the  input  to  the  feature-based 
analysis  code. 


Due  to  the  rather  broad  set  of  possible  interpretations  of  the  term  “feature”,  each  kind  of  topological  feature 
extraction  is  itself  considered  an  application  program  and,  therefore,  is  not  a  standard  part  of  the  interrogation 
library.  The  process  of  topological  feature  extraction  is  currently  programmed  using  the  ray-tracing  interrogation 
features  described  above.* 


3.3  3-D  Surface  Mesh 

Combinatorial  Solid  Geometry  (CSG)  models  are  formed  by  the  boolean  combination  of  “primitive”  solids.7 
For  example,  a  plate  with  a  hole  is  most  easily  modeled  as  a  plate  primitive  minus  a  cylinder  primitive.  It  is 
important  that  in  CSG  models,  there  is  no  explicit  representation  of  the  surfaces  of  the  solids  stored;  indeed,  for 
complex  boolean  combinations  of  complex  primitives,  some  of  the  resultant  shapes  may  have  very  convoluted 
topology  and  surfaces  that  may  be  at  best  high  degree  polynomials. 

There  are  many  applications  that  would  benefit  from  being  able  to  express  an  approximation  of  these  complex 
shapes  created  using  CSG  modeling  as  a  collection  of  planar  N-sided  polygons  (N-gons)  which  together  enclose 
roughly  the  same  volume  of  space  as  the  original  CSG  solid.  The  most  obvious  such  application  is  to  drive 
polygon-based  rendering  routines  (lighting  modules)  for  predictive  optical  signatures.  On  many  modern 
workstations  there  is  direct  hardware  or  firmware  support  for  high-speed  rendering  of  polygons.  In  addition,  there 
are  whole  collections  of  polygon-based  predictive  infrared  and  radar  signature  programs.  Some  of  the  most 
accurate  radar  signatures  have  been  calculated  using  the  Method  of  Moments.17  This  approach  requires  a  three 
dimensional  surface  tessellation  to  sub-wa  'elength  resolution  over  an  entire  geometric  model. 

A  sensible  strategy  for  converting  a  CSG  model  to  the  equivalent  approximate  three  dimensional  surface  mesh 
is  to  tackle  the  problem  in  two  parts.  First,  a  routine  has  to  be  written  to  convert  each  of  the  primitive  solids  into 
tessellated  form.  Second,  a  routine  has  to  be  written  to  take  two  tessellated  objects  and  combine  them  according 
to  a  boolean  operation  (union,  intersection  or  subtraction)  back  into  a  consistent  set  of  solid  tessellated  objects. 
Until  very  recently,  this  second  step  has  proven  extremely  difficult.  The  topology  of  solid  tessellated  objects  has 
traditionally  been  represented  using  the  “winged-edge”  data  structure.  Within  the  winged-edge  representation,  an 
edge  represents  the  boundary  or  intersection  between  exactly  two  faces.  Unfortunately,  this  structure  fails  to 
handle  other  valid  configurations,  such  as  an  edge  being  shared  by  four  faces,  or  an  edge  being  part  of  only  a  single 
(“dangling”)  face.  These  other  configurations  arise  when  the  topology  of  an  object  is  not  that  of  a  simple  3- 
manifold,  i.e.  when  the  topology  of  the  object  can  not  be  mapped  to  a  sphere.  These  non-3-manifold  conditions 
arise  in  the  construction  of  finite  element  meshes,  and  from  the  use  of  the  boolean  intersection  operation. 

The  inability  of  the  winged-edge  data  structure  to  represent  non-3-manifold  configurations  resulted  in  the 
development  of  a  more  general  data  structure18,19.  This  new  data  structure  has  been  dubbed  alternately  the 
“radial-edge”,  “Non-Manifold  Geometry”  (NMG),  or  “n-Manifold  Geometry”  data  structure.  The  radial-edge 
representation  provides  topological  links  between  all  faces  which  share  an  edge.  This  single  representation  has  the 
ability  to  handle  n-manifolds  (Mn)  for  0  <  n  <  3:  3-manifolds  (solid  objects),  2-manifolds  (lone  faces,  not  part  oFa 
solid),  1-manifolds  (lone  edges,  not  part  of:  a  face)  and  O-manifolds  (lone  points,  not  part  of  an  edge). 


t  See  Ref.  5,  Figs.  8  and  0. 


10 


Consider  the  intersection  of  two  solids  that  share  only  a  single  face,  edge,  or  vertex.  The  result  of  the 
intersection  will  be  a  manifold  object  of  dimension  less  than  three.  The  winged-edge  data  structure  is  unable  to 
represent  these  boolean  results,  because  it  requires  every  edge  to  adjoin  exactly  two  faces.  Thus,  the  winged-edge 
representation  is  not  closed  under  boolean  operations.  To  overcome  this  lack  of  closure,  winged-edge  systems 
substitute  regularized  boolean  operators  which  are  defined  to  produce  only  3-manifold  results;  all  lesser  dimension 
results  are  discarded.  In  contrast,  because  all  manifold  objects  of  dimension  less  than  three  can  be  represented 
using  the  NMG  data  structures,  NMG  objects  are  closed  under  boolean  operations. 

Employing  the  NMG  representation  for  faceted  solid  objects  gives  rise  to  the  rich  set  of  possibilities  diagramed 
in  Fig.  4.  From  this  diagram  it  should  be  clear  that  the  final  evaluated  NMG  solid  object  can  be  employed  in  a 
variety  of  ways.  The  primary  use  will  be  for  input  to  analysis  codes  that  need  an  approximate  3-D  surface  mesh  of 
the  solid  model.  In  this  case,  the  NMG  objects  are  sent  across  the  interface,  either  directly  into  an  application,  or 
via  a  triangulator  that  turns  the  planar  N-gon  faces  of  the  NMG  objects  into  simple  triangle  lists  and  thence  to 
applications  such  as  SARSIM.  However,  a  very  powerful  second  use  will  be  to  create  new  faceted  shapes  which  are 
then  stored  back  in  the  database  as  new  geometric  objects,  suitable  for  future  editing  or  analysis. 

While  a  detailed  description  of  the  NMG  data  structures  is  beyond  the  scope  of  this  paper,  there  are  several 
advantageous  properties  of  the  NMGs  that  are  worth  mentioning.  The  NMG  representation  maintains  full 
topology  information,  so  that  the  relationships  between  vertices,  edges,  loops,  faces  and  shells  are  continuously 
available.  The  geometry  information  associated  with  a  planar  face  is  the  plane  equation  (which  includes  the 
outward-pointing  surface  normal);  the  plane  equation  does  not  have  to  be  re-derived  from  the  vertices.  For 
applications  that  would  prefer  visual  realism  rather  than  geometric  fidelity,  there  is  room  in  the  vertex  geometry 
structure  to  carry  around  a  “phony”  normal  for  each  vertex,  suitable  for  intensity  interpolation  shading 
algorithms  used  in  Gouraud  shading,"1  or  for  normal-vector  interpolation  shading  algorithms. 

One  of  the  most  exciting  current  research  projects  at  BRL  is  the  extension  of  the  NMG  framework  to  permit 
faces  either  to  be  planar  N-gons,  or  trimmed  Non-Uniform  Rational  B-Splines  (“trimmed  NURBS”).  This  will 
permit  many  of  the  tessellation  operations  to  be  implemented  exactly,  rather  than  as  approximations.  This  will 
also  permit  solids  to  enjoy  the  economy  of  having  most  faces  be  represented  as  planar  N-gons,  which  are  very 
compact  and  efficient  to  process,  while  those  few  faces  that  require  sculptured  surface  shape  control  can  be 
represented  as  trimmed  NURBS.  This  combination  provides  both  efficiency  and  full  shape  control  in  the  rich  n- 
manifold  topological  framework;  a  combination  that  does  not  exist  in  any  current  commercial  CAD  system. 

3.4  3-D  Volume  Mesh 

Many  forms  of  energy  flow  analysis,  such  as  heat  flow,  vibrational  analysis  (acoustic  energy  flow),  and  stress 
analysis  require  the  use  of  3-D  Finite-Element  Mesh  (FEM)  techniques.  While  there  has  been  some  work  on  using 
the  ray-tracing  paradigm  to  construct  finite  element  and  finite  difference  meshes22  it  has  been  difficult  to  deal  with 
high  spatial  frequency  (fine  detail)  portions  of  the  model.  In  particular,  meshing  small  diameter  pipes  is 
problematic;  undersampling  can  cause  the  pipe  to  incorrectly  be  separated  into  multiple  pieces. 

In  order  to  improve  on  the  current  state  of  affairs,  it  seems  necessary  to  provide  support  for  the  generation  of 
volume  meshes  directly  as  part  of  the  application  interface.  This  would  provide  the  meshing  algorithm  to  have 
unrestricted  access  to  the  underlying  geometry,  the  space  partitioning  tree,  and  other  internal  data  in  order  to 
perform  a  better  job. 

Even  more  promising  still  would  be  a  strategy  that  takes  advantage  of  the  NMG  support.  A  first  pass  might 
tessellate  the  model  and  evaluate  the  booleans  to  produce  a  surface  mesh.  The  second  pass  would  then  take  the 
surface  mesh  and  fill  the  interior  (or  exterior)  volumes  with  appropriately  chosen  volume  elements.  A  very  good  fit 
could  probably  be  achieved  using  only  parallelepiped  (“brick”)  elements  and  20-node  “superelements”.  The  brick 
elements  would  be  used  to  fill  interior  volume  that  does  not  border  on  a  face,  and  the  superelements  would  be  used 
for  volume  that  contacts  a  face.  Recourse  could  be  made  back  to  the  underlying  geometry  (perhaps  via  firing  a  few 
well  chosen  rays)  to  get  the  curvature  of  the  superelement  faces  to  match  the  curvature  of  the  underlying 
primitive,  rather  than  having  to  rely  strictly  on  the  NMG  planar-face  approximation. 
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Figure  4.  Logic  flow  for  n-Manifotd  Geometry  (NMG)  processing. 


3.5  3-D  Volume  Elements  (Voxels) 

A  representation  ’.vhich  is  similar  to  the  finite-element  mesh  is  based  on  Volume  Elements  (VOXELS).  There 
are  two  distinct  kinds  of  voxels.  The  first  kind  of  voxels  can  be  considered  a  special  case  of  volume  meshing 
discussed  previously,  in  which  the  model  is  “diced”  into  a  large  collection  of  homogeneous  parallelepiped  (“brick”) 
elements.  As  one  example,  ERIM  has  a  utility  program  which  uses  ray-tracing  to  convert  BRL  CSG-format 
geometry  to  this  kind  of  voxel  representation  to  feed  a  first-principles  IR  model.23  In  contrast  to  the  ERIM’s  voxel 
modeling  approach,  the  Physically  Reasonable  Infra-Red  Signature  Model  (PRISM),  predictive  IR  code,  developed 
at  the  Keweenaw  Research  Center,  is  based  on  a  geometry  of  flat  plates.24  Two  codes  linking  the  BRL  CSG-format 
geometry  to  PRISM  have  been  developed:  the  TACOM  FRED  editor”5  and  the  BRL  irprep  program.26 

A  distinctly  different  form  of  voxel  representation  is  based  upon  the  use  of  8-way  binary  space  subdivision 
stored  using  an  “oct-tree”  data  structure.  In  this  technique,  the  model  is  enclosed  in  a  bounding  box.  The 
bounding  box  is  evenly  split  along  the  X,  Y  and  Z  axes  to  form  eight  smaller  boxes.  This  algorithm  is  applied 
recursively  so  that  all  boxes  which  are  neither  entirely  full  nor  entirely  empty  are  repeatedly  split,  until  the  size  of 
the  voxels  satisfies  some  termination  condition.  In  this  way,  small  voxels  that  lie  along  the  surface  of  objects  can 
fit  arbitrarily  tightly  to  the  surface,  while  the  interior  of  an  isomorphic  region  will  be  contained  primarily  in  a 
single  large  voxel. 

The  oct-tree  representation  provides  the  application  program  with  a  homogeneous  geometric  representation 
based  entirely  on  cubes  of  varying  size.  Having  such  a  homogeneous  representation  can  often  greatly  ease  the  task 
of  algorithm  development.  On  the  other  hand,  achieving  a  good  approximation  of  curved  objects  using  cubes 
requires  a  huge  number  of  voxels  to  be  used,  resulting  in  very  large  voxel  datasets,  and  an  exponential  increase 
(order  N8)  in  the  number  of  element-to-element  equations  to  be  solved.  The  oct-tree  approach  to  IR  signature 
generation  is  employed  in  the  BRL-CAD  program  Igt."7 
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3.0  Homogeneous  Trimmed  B-Splines 

When  support  for  trimmed  NURB  faces  has  been  added  to  the  NMG  capability,  it  will  be  possible  to  represent 
all  existing  primitives  either  with  exact  rational  B-spline  versions,  or.  with  very  good  rational  B-spline 
approximations.  This  could  be  done  even  for 'faces  that  were  completely  planar. 

This  offers  the  hope  that  it  might  be  possible  (albeit  memory  intensive)  to  convert  an  entire  CSG  solid  model 
into  a  homogeneous  collection  of  non-uniform  rational  B-spline  faces  organized  in  a  n-manifold  topological  data 
structure.  In  addition  to  the  conceptual  simplicity  afforded  by  having  a  uniform  representation  for  shape,  this 
affords  the  opportunity  to  create  new  analysis  codes  that  can  process  curved  surfaces,  yet  at  least  initially  only 
have  to  deal  with  one  kind  of  shape.  This  would  also  provide  a  very  direct  and  natural  interface  to  spline-based 
and  Bezier-patch28  based  modeling  systems. 

3.7  Analytic  Analysis 

Given  a  homogeneous  geometric  representation  such  as  the  Trimmed  B-Splines  just  discussed  which  also  has  an 
analytic  representation,  a  further  processing  capability  arises.  Rather  than  interrogating  the  data  base  by  means 
of  sampling  or  subdivision  techniques,  the  direct  mathematical  manipulation  of  the  source  geometry  through  its 
parametric  representation  becomes  possible.  Calculations  of  physical  properties  requiring  integration  over  a  surface 
can  often  be  evaluated  with  greater  accuracy  using  an  explicit  analytic  calculation  than  could  be  provided  by 
numerical  methods.  While  this  may  be  difficult  in  general  due  to  the  complexity  of  a  parametric  expression,  some 
classes  of  surface  representations  may  be  good  candidates  for  this  approach.  Splines,  for  example,  arc  piecewise- 
polynomial  functions  which  have  relatively  simple  Fourier  transform  representations.  Since  2-D  spatial  Fourier 
transforms  arise  frequently  in  far-field  electromagnetic  scattering  calculations,  exploitation  of  the  parametric  spline 
form  is  of  interest  in  predictive  scattering  calculations.  Direct  use  of  spline  parameters  in  a  PO  scattering  model  is 
part  of  the  methodology  used  at  the  Aircraft  Division  of  The  Northrop  Corporation.29 

Support  for  this  tack  may  be  inferred  from  the  work  of  Schneider  and  Peden  who  have  exploited  analytic 
methods  for  calculation  of  radar  cross  sections  using  dielectric  ellipsoids30  in  the  detection  of  buried  targets.31  Here 
the  approach  involves  the  approximation  of  the  geometry  of  interest  with  a  shape  for  which  the  analytic  solution  of 
the  electromagnetic  scattering  problem  can  be  solved.  Great  care  must  be  taken  to  insure  that  the  errors 
introduced  by  the  geometry  are  a  small  perturbation  of  the  desired  solution. 

With  the  rapidly  developing  pote  lal  of  symbolic  calculation,  treatment  of  seemingly  impossible  formulas 
resulting  from  the  geometry/physics  ir‘eraction  may  become  tenable.  This  could  help  to  reduce  the  trend  towards 
employing  numerical  methods  at  the  onset  of  a  problem  and  avoid  some  of  the  accompanying  instabilities  and 
errors. 


4.  SUMMARY 

In  this  paper  we  have  discussed  some  issues  basic  to  the  prediction  of  signatures  in  support  of  Automatic  Object 
Recognition,  in  particular  the  way  in  which  three-dimensional  geometry  and  material  data  are  linked  to  certain 
applications.  This  point  was  illustrated  via  a  discussion  of  three  predictive  radar  codes,  each  designed  for 
essentially  the  same  application,  but  each  nevertheless  employing  unique  geometric  methods.  The  descriptive 
approach  and  the  manner  in  which  each  is  linked  to  the  physics  of  the  codes »has  been  shown  to  both  enable  and 
constrain  algorithmic  exploitation. 

The  methods  for  linking  geometry  to  applications  codes  were  described  in  turn.  The  BRL-CAD  package  was 
used  to  illustrate  each  approach.  As  described  in  the  APPENDIX,  the  strategy  utilized  in  BRL-CAD  is  based  on 
an  inhomogeneous  collection  of  closed-surface  geometries,  of  which  the  variety  of  shapes  is  constantly  expanding. 
By  this  tack,  in  addition  to  using  geometry  files  created  with  BRL-CAD,  three-dimensional  geometry  files  that  were 
developed  under  other  systems  can  either  be  mapped  into  the  BRL-CAD  data  base  exactly  (if  the  corresponding 
geometric  forms  already  exist)  or  the  data  base  can  be  extended  to  support  any  important  new  forms  (so  that  the 
conversion  becomes  an  exact  mapping  as  well).  By  this  method  the  enormous:  costs  of  geometry  generation  can  be 
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recovered;  any  model  developed  on  a  true  3-D  modeling  system  can  be  transferred  through  file  importation  to  the 
general  and  open  environment  of  BRL-CAD.* 

Another  important  point  that  was  emphasized  was  the  possible  utility  of  each  particular 'interface.  A  related 
goal,  therefore,  of  BRL-CAD  has  been  to  support  a//  such  interfaces  from  the  extensible  data  base  of  geometry 
types.  Many  of  the  forms  have  been  implemented;  the  NMGs  which  will  provide  homogeneous  polygon  support  ore 
in  the  final  stages  of  development. 

By  these  developments  it  is  our  intent  to  support  the  AOR  community  to  the  maximum  extent  possible  with 
the  tools  needed  to  develop  the  multi-spectral,  multi-sensor  simulations  critical  to  autonomous-sensor  applications. 

APPENDIX 

The  solid  modeling  package  called  BRL-OAD9’15  was  originally  developed  to  support  the  input  requirements  of 
vulnerability  and  lethality  modeling  at  the  Ballistic  Research  Laboratory.  Its  capabilities  are  briefly  summarized 
here: 

•  BRL-CAD  is  composed  of  more  than  200,000  lines  of  portable  C  language  source  code  which  support: 

□  Solid  geometric  editor  (mged) 

□  Ray  tracing  utilities 

□  Lighting  models  for  a  variety  of  frequencies 

□  Many  image-handling,  data-comparison,  image-processing  and  other  supporting  utilities 

•  The  set  of  closed-surface  (inhomogeneous)  geometrical  ^presentations  supported  by  BRL-CAD  include: 

□  The  original  Constructive  Solid  Geometry  (CSG)  BRL  data  base 

□  Non-Uniform  Rational  B-Spline  Surfaces  (NURBsj 

□  The  faceted  data  representation  (PATCH)  developed  by  Falcon/  Denver  Research  Institute  and  used  by  the 
Navy  and  Air  Force  for  vulnerability  calculations. 

•  It  supports  association  of  material  (and  other  attribute  properties)  with  geometry  which  is  critical  to 
subsequent  applications  codes. 

•  It  supports  a  set  of  extensible  interfaces  by  means  of  which  geometry  (and  attribute  data)  are  passed  to 
applications: 

□  Ray  casting 

□  Topological  representation 

□  3-D  Surface  Mesh  Generation 

□  3-D  Volume  Mesh  Generation 

□  Analytic  (Homogeneous  Spline)  representation 

•  Source  code  for  BRL-CAD  has  been  distributed  to  more  than  650  computer  sites  world  wide  throughout 
Government,  Industry  and  Academia. 

•  In  addition  to  the  vulnerability  and  signature  codes  generated  by  the  BRL,  many  BRL-CAD  based  applications 
codes  have  been  built  by  others,  including  applications  developed  by  workers  at  TACOM/  Keweenaw  Research 
Center,  ERIM,  Northrop,  MITRE,  University  of  Illinois  and  scores  of  other  sites. 

Figure  A1  gives  a  general  layout  of  the  BRL-CAD  database  representation.  A  number  of  fundamental 
strategies  are  key  here.  First,  the  primary  data  base  is  inhomogeneous;  that  is,  many  types  of  geometric 
representations. are  allowed  so  long  as  they  represent  fully  enclosed  space.  New  shapes  are  added  to  the  geometric 


t  An  important  caveat  here  is  that  the  mathematical  form  of-the;geometric  representation  and  data  base  to  be  imported  must  be  known. 
Unfortunately,  most  commercial  CAD  vendors  are  unwilling  to  provide  this  information  for  files  generated  by  their  systems.  since  that 
might  free  the  user  from  "vendor  lockin". 
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Figure  Al.  Org&niiation  of  BRL-CAD  Database  Representation. 


Figure  A2.  Current  Applications  of  BRL-CAD.1*^’10’16 
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data  base  from  time  to  time  and  the  BRL-CAD  package  has  a  number  of  utilities  written  expressly  to  aid  that 
process. 

Second,  sppcific  interface  standards  have  been  adopted  at  the  points  where  geometry  is  normally  passed  to 
applications  ( e.g .  ray  casting,  topology,  surface  mesh).  By  this  strategy  changes  which  take  place  in  the  geometric 
data  base  itself  are  isolated  from  the  specific  means  of  geometry  linkage.  By  means  of  this  strategy  both  (a)  an 
extremely  large  primary  data  base  of  geometry  types  can  be  supported  in  a  data-storage  efficient  fashion  and  (b) 
various  application  codes  can  link  to  that  data  base  in  the  manner  most  suitable  for  the  application. 

As  noted  above,  although  the  original  application  of  BRL-CAD  was  in  support  of  vulnerability  and  lethality 
analyses,  its  utility  has  been  extended  much  beyond  those  applications.1’7,10,  5  A  graphical  depiction  of  current 
uses  is  given  in  Fig.  A2. 
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ATTN:  AMSEL-RD-NV-V  (John  llo) 

Fort  Belvoir,  VA  22060-5677 

1  Director 

Center  for  Night  Vision  and  Elcctro-Optics 
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Warren,  MI  48397-5000 

1  Commander 

US  Army  Tank-Automotive  Command 
ATTN:  AMSTA-CV  (COL  Becking) 
Warren,  MI  48397-5000 

2  Commander 

US  Army  Tank-Automotive  Command 
ATTN:  AMSTA-NKS  (D.  Cyaye) 

(J.  Rowe) 

Warren,  MI  48397-5000 

2  Commander 

US  Army  Tank-Automotive  Command 
ATTN:  AMSTA-RG  (R.  Munt) 

(R.  McClelland) 

Warren,  MI  48397-5000 


3  Commander 

US  Army  Tank-Automotive  Command 
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US  Army  Tank-Automotive  Command 
ATTN:  AMSTA-RSIC  (Sam  Goodman) 
Warren,  MI  48090-5000 

1  Commander 

US  Army  Tank-Automotive  Command 
ATTN:  AMSTA-VS  (Brian  Bonkosky) 
W arren,  MI  48090-5000 

6  Commander 

US  Army  Tank-Automotive  Command 
ATTN:  AMSTA-ZE  (R.  Asoklis) 

AMSTA-ZEA  (C.  Robinson) 

(R.  Gonzalez) 
AMSTA-ZS  (D.  Rees) 
AMSTA-ZSS  (J.  Thompson) 

(J.  Soltez) 

Warren,  MI  48397-5000 

1  Commander 
HQ,  TRADOC 

ATTN:  Asst  Dep  Chief  of  Staff 
for  Combat  Operations 
Fort  Monroe,  VA  23651-5000 

2  Director 

HQ,  TRAC  RPD 

ATTN:  ATRC-RP  (COL  Brinkley.) 

ATRC-RPR  (Mark  W.  Murray.) 
Ft.  Monroe,  VA  23651-5143 

1  Director 

US  Army  Cold  Regions  Research  and 
Development  Laboratory 
ATTN:  Technical  Director  (Lewis  Link) 
72  Lyme  Road 
Hanover,  NH  03755 


24 


No.  of 

Copies  Qrga.ni7a.ti0n 


Organization 


No.  of 
Copies 


l  US  Army  Corps  of  Engineers 

Assistant  Director  Research  and  Development 
Directorate 
ATTN:  Mr.  B.  Benn 
20  Massachusetts  Avenue,  NW 
Washington,  DC  20314-1000 

1  Commander 

US  Army  Operational  Test  and  Evaluation 
Agency 

ATTN:  MG  Stephenson 
4501  Ford  Avenue 
Alexandria,  VA  22302-1458 

1  Commander 

US  Army  Vulnerability  Assessment 
Laboratory 

ATTN:  SLCVA-CF  (Gil  Apodaca) 

White  Sands  Missile  Range,  NM  88002-5513 

1  Director 

TRAC-WSMR 

ATTN:  ATRC-RD  (McCoy) 

WSMR,  NM  88002-5502 

2  US  General  Accounting  Office 
Program  Evaluation  and  Methodology 

Division 

ATTN:  Robert  G.  Orwin 
Joseph  Sonnefeld 
Room  5844 
441  G  Street,  NW 
Washington,  DC  20548 

1  Director 

Office  of  the  Deputy  Under  Secretary 
of  the  Army,  Operations  Research 
Study  Program  Management  Agency 
ATTN:  SFUS-SPM/E.  Visco 
Washington,  DC  20310-0102 

1  Director 

US  Army  Industrial  Base  Engineering  Activity 
ATTN:  AMXIB-MT 
Rock  Island,  IL  61299-7260 


1  Director 

US  Army  Industrial  Base  Engineering  Activity 
ATTN:  AMXIB-PS  (Steve-McGlone) 

Rock  Island,  IL  61299-7260 

3  Director 

US  Army  Engineer  Waterways  Experiment 
Station 

ATTN:  WESEN  (Dr.  V.  LaGarde) 

(Mr.  W.  Grabau) 

WESEN-C  (Mr.  David  Meeker) 

PO  Box  631 

Vicksburg,  MS  39180-0631 

1  US  Army  Engineer  Topographic  Laboratories 
ATTN:  Technical  Director  (W.  Boge) 

Fort  Belvoir,  VA  22060-5546 

1  Commander 

US  Army  Operational  Test  and  Evaluation 
Agency 

ATTN:  LTC  Gordon  Crupper 
4501'Ford  Ave.  #870 
Alexandria,  VA  22302-1435 

1  Lawrence  Livermore  National  Laboratories 
PO  Box  808  (L-3321) 

ATTN:  MarkV/ilkins 
Livermore,  CA  94551 

3  Los  Alamos  National  Laboratories 
ATTN:  MS  985,  Dean  C.  Nelson 
MS  F600,  Gary  Tietgen 
MS  G787,  Terrence  Phillips 
PO  Box  1663 
Los.Alamos,  NM  87545 

1  Los  Alamos  National  Laboratories 

ATTN:  MS  F681,  LTC  Michael  W.  Ziehmn 

USMC 

POBox  1668 

Los  Alamos,  NM  87545 

1  Sandia  National  Laboratories 
Department  913 
ATTN:  Ron  Andreas 
Albuquerque,  NM  87185-5800 
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1  Sandia  National  Laboratories 
Division  1611 
ATTN:  Tom  James 
Albuquerque*  NM  87185 

1  Sandia  National  Laboratories 
Division  1623 
ATTN:  Larry  Hostetler 
Albuquerque,  NM  87185 

1  Sandia  National  Laboratories 
ATTN:  Gary  W.  Richter 
PO  Box  969 
Livermore,  CA  94550 

1  Commander 

US  Naval  Air  Systems  Command 
JTCG/AS  Central  Office 
ATTN:  5164J  (LTC  James  B.  Sebolka) 
Washington,  DC  20361 

1  Commander 

ADR  Program  Manager 
CODE  AIR-411121 
ATTN:  Tom  Furlough 
Naval  Air  Systems  Command 
Washington,  DC  20361-4110 

1  Commander 

US  Naval  Ocean  Systems  Center 
ATTN:  Earle  G.  Schweizer 
Code  000 

San  Diego,  CA  92151-5000 
4  Commander 

US  Naval  Surface  Warfare  Center 
ATTN:  Gregory  J.  Budd 
James  Ellis 
Barbara  J.  Harris 
Constance  P.  Rollins 
Code  G13 

Dahlgren,  VA  22448-5000 
1  Commander 

US  Naval  Surface  Warfare  Center 
ATTN:  Frank  Fassnacht 
10901  New  Hampshire  Ave. 

Code  N15 

Silver  Spring,  MD  20903-5000 


1  Commander 

US  Naval  Surface  Warfare  Center 
ATTN:  Norma  D.  Holland 
Code  R-14 

Silver  Spring,  Md  20903-5000 
1  Commander 

US  Naval  Surface  Warfare  Center 
ATTN:  Dr.  F.E.  Baker 
10901  New  Hampshire  Ave. 

Silver  Spring,  MD  20903-5000 

1  Commander 

US  Naval  Surface  Warfare  Center 
ATTN:  William  Embcrson 
Code  H021 

10901  New  Hampshire  Ave. 

Silver  Spring,  MD  20903-5000 

1  Commander 

US  Naval  Surface  Warfare  Center 
ATTN:  M.  John  Timo 
10509  Edgefield  Drive 
Adelphi,  MD  20783-1130 

2  Commander 

US  Naval  Weapons  Center 
ATTN:  Ed  Patterson 

Dr.  Helen  Wang 
Code  3313 

Bldg  1400,  Room  B17 
China  Lake,  CA  93555 

1  Commander 

US  Naval  Weapons  Center 
ATTN:  David  H.  Hall 
Code  3181 

China  Lake,  CA  93555-6001 

1  Commander 

US  Naval  Weapons  Center 
ATTN:  Mark  D.  Alexander 
Code  3894 

China  Lake,  CA  93556-6001 
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1  Commander 

US  Naval  Weapons  Center 
ATTN:  Melvin  H.  Keith 
Code  39104 

China  Lake,  CA  93555 

1  Commander 

US  Naval  Weapons  Center 
ATTN:  Tim  Horton 
Code  3386 

China  Lake,  CA  93555 

1  Commander 

US  Naval  Weapons  Center 
ATTN:  Robert  Cox 
Code  3917 

China  Lake,  CA  93555-6001 

1  Commander 

US  Naval  Civil  Eng  Laboratories 
ATTN:  John  M.  Ferritto 
Code  L53 

Port  Hueneme,  CA  93043 

1  Naval  Postgraduate  School 
ATTN:  Dr.  Robert  E.  Ball 
642  Toyon  Drive 
Monterey,  CA  93940 

1  Naval  Postgraduate  School 
ATTN:  Dr.  Michael  J.  Zyda 
Department  of  Computer  Science 
Code  52 

Monterey,  CA  93943-5000 

1  Naval  Postgraduate  School 

Department  of  National  Security 
ATTN:  Dr.  Joseph  Sternberg 
Code  73 

Monterey,  CA  93943 

1  Commander 

US  Naval  Sea  Systems  Command 
ATTN:  William  G.  Boyce 
Code  56Y52 
Washington,  Dc  20362 


No.  of 

Copies  Organization 

1  Commander 

US  Naval  Sea  Systems  Command 
ATTN:  Granville  W.  Broome 
SEA  5011 

2521- Jefferson  Davis  Hwy. 
Arlington,  VA  22202 

1  Commander 

US  Naval  Sea  Systems  Command 
ATTN:  Philip  M.  Covich 
SEA  55X 

Washington,  DC  20362-5101 

2  Commander 

US  Naval  Sea  Systems  Command 
ATTN:  CPT  Charles  Calvano  USN 
Robert  Keane,  Jr. 

SEA  50 

Washington,  DC  20362-5101 
1  Commander 

US  NavallSea  SystemsiCommand 
ATTN:  Oliver  F.  Braxton 
2521  Jefferson  Davis  Hwy. 
Arlington,  VA  22202 

1  Commander 

US  Naval  Sea  Systems  Command 
ATTN:  Donald  Ewing 
Code  503 

2521  Jefferson  Davis  Hwy. 
Arlington, 'VA  22202 

1  Commander 

US  Naval  Sea  Systems  Command 
ATTN:  Larrie  D.  Ferreiro 
SEA  501 

2521  Jefferson  Davis  Hwy. 
Arlington,  VA  22202 

1  Commander 

US  Naval  Sea  Systems  Command 
ATTN:  Anthony  F.  Johnson 
SEA  05R2 

Washington,  DC  20362-5101 
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1  Commander 

US  Naval  Sea  Systems  Command 
ATTN:  CPT  William  E.  Mahew  USN 
PMS  423 

Washington,  DC  20362-5101 

1  Commander 

US  Naval  Sea  Systems  Command 
ATTN:  Carl  H.  Pohler 
Code  05R23 

Washington,  DC  20362-5101 

1  Commander 

US  Naval  Sea  Systems  Command 
ATTN:  Ronald  P.  Kramer 
SEA  50143 

2521  JefTerson  Davis  Hwy. 

Arlington,  VA  22202 

1  Commander 

US  Naval  Sea  Systems  Command 
ATTN:  CPT  R.  Percival  USN 
SEA05T 

2521  Jefferson  Davis  Hwy. 

Arlington,  VA  22202 

1  Commander 

US  Space  and  Naval  Warfare  Systems 
Command 

ATTN:  Paul  Wessel 
Code  30T 

Washington,  DC  20363-5100 
1  Commander 

Intelligence  Threat  Analysis  Center 
ATTN:  PSD-GAS/John  Bickle 
Washington  Navy  Yard 
Washington,  DC  20374 

1  Commander 

Intelligence  Threat  Analysis  Center 
ATTN:  Ron  Demeter 
Washington  Navy  Yard,  B-213,  Stop  314 
Washington,  DC  20374 


No.  of 

Copies  Organization 

1  Commander 

Intelligence  Threat  Analysis  Center 
ATTN:  Tim  Finnegan 
Washington  Navy  Yard,  B-213 
Washington,  DC  20374 

2  Commander 

Intelligence  Threat  Analysis  Center 
Intell  Image  Prod  Div 
ATTN:  John  Creighton 
A1  Fuerst 

Washington  Navy  Yard,  Bldg  213  (IAX-0-II) 
Washington,  DC  20374 

2  Commander 

David  W.  Taylor  Naval  Ship  and 
Development  Center 
ATTN:  W.  Conley 
J.  Schot 

Bethesda,  MD  20084 

1  Commander 

Eglin  Air  Force  Base 

AD/ENL 

ATTN:  Robert  L.  Stovall 
Eglin  AFB,  FL  32542 

1  Commander 

USAF  HQ  ESD/PLEA 
,  Chief,  Engineering  and  Test  Division 
ATTN:  Paul  T.  Courtoglous 
Hanscom  AFB,  MA  01730 

2  Commander 

AFATL 

ATTN:  AGA  (Lawrence  Jones) 

(Mickie  Phipps) 

Eglin  AFB,  FL  32542-5434 

1  Commander 

AFEWG 

ATTN:  AFEWC/SAXE  (Bod  Eddy) 

Keljy  AFB,  TX  78243-5000 

1  Commander 
AFVv’AL  /'AARA 
ATTN:  Ed  Zelano 
Wright-Patterson  AFB,  OH  45433 
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1  Commander 
AFWAL/FIES 
ATTN:  James  Hodges  Sr. 
Wright-Patterson  AFB,  OH 
45433-6523 


2  Commander 
AFWAL/MLTC 
ATTN:  LT  Robert  Carringer 
Dave  Judson 

Wright-Patterson  AFB,  OH 
45433-6533 


1  Commander 
ASB/XRM 

ATTN:  Gerald  Bennett 
Martin  Lentz 

Wright-Patterson  AFB,  OH 

1  Commander 
WRDC/AARA 
ATTN:  Michael  L.  Bryant 
Wright-Patterson  AFB,  OH  45433 

1  Commander 
FTD/SDMBA 
ATTN:  Charles  Darnell 
Wright-Patterson  AFB,  OH  45433 

1  Commander 
FTD/SDMBU 
ATTN:  Kevin  Nelson 
Wright-Patterson  AFB,  OH  45433 

1  Commander 
FTD/SQDRA 
ATTN:  Greg  Koesters 
Wright-Patterson  AFB,  OH 
45433-6508 

1  Commander 
FTD 

ATTN:  Tom  Reinhardt 
Wright-Patterson  AFB,  OH  45433 

1  -Commander 
FTD/SCRS 

ATTN:  Amy  Fox  Schalle 
Wright-Patterson  AFB,  OH  45433 


1  Commander 
FTD/SDJEO 
ATTN:  Robert  Schalle 
Wright-Patterson  AFB,  OH  45433 

1  Commander 
FTD/SDAEA 
ATTN:  Joe  Sugrue 
Wright-Patterson  AFB,  OH  45433 

1  Commander 
AFWAL/AARA 
ATTN:  Vincent  Velten 
Wright-Patterson  AFB,  OH  45433 

1  Commander 
FTD/SQDRA 
ATTN:  Larry  E.  Wright 
Wright-Patterson  AFB,  OH  45433 

1  Commander 
AD/CZL 

ATTN:  James  M.  Heard 
Eglin  AFB,  FL  32542-5000 

1  Commander 
AD/ENY 

ATTN:  Dr.  Stewart  W.  Turner 
Director  of  Engineering  Analysis 
Eglin  AFB,  FL  32542-5000 

2  Commander 
AD/ENYW 

ATTN:  2LT  Michael  Ferguson 
Jim  Richardson 
Eglin  AFB,  FL  32542-5000 

1  Commander 

Air  Force  Armament  Laboratory 
ATTN:  AFATL/DLY  (James  B.  Flint) 
Eglin  AFB,  FL  32542-5000 
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4  Commander 
US  Army  FSTC 
ATTN:  Greg  Crawford 
David  P.  Lutz 
Suzanne  Hall 
Charles  Hutson 
220  Seventh  Avenue 
Charlosttesville,  YA  22901-5396 


1  Commander 

US  Army  FSTC/CA3 
ATTN:  Scott  MingledorlT 
220  Seventh  Avenue 
Charlottesville,  VA  22901-5396 


1  Commander 

US  Army  FSTC  (UK) 

ATTN:  MAJ  Nigel  Williams 
220  Seventh  Avenue 
Charlottesville,  VA  22901-5396 


1  Commander 
US  Army  FSTC 
ATTN:  Dr.  Tim  Small 
220  Seventh  Avenue 
Charlottsville,  VA  22901-5396 


1  Defense  Intelligence  Agency 
ATTN:  DB-6E3  (Jay  Hagler) 
Washington,  DC  20340-6763 


6  Institute  for  Defense  Analyses  (IDA) 
ATTN:  Mr.  Irwin  A.  Kaufman 
Mr.  Arthur  0.  Kresse 
Mr.  Arthur  Stein 
Dr.  Lowell  Torinessen 
Mr.  Benjamin  W.  Turner 
Ms.  Sylvia  L.  Waller 
1801  N.  Beauregard  Street 
Alexandria,  VA  22311 


1  Institute  for  Defense  Analyses 
ATTN:  Carl  F.  Kossaclc 
1005  Athens  Way 
Sun  City,  FL  33570 


No.  of 

Copies  Organization 

1  Institute  for  Defense  Analyses 

ATTN:  Dr.  Natarajan  Subramonian 
14309  Hollyhock  Way 
Burtonsville,  MD  20866 

1  Department  of  Commerce 

National  Institute  of  Standards  and 
Technology 

Manufacturing  Systems  Group 
ATTN:  B.  Smith 
Washington,  DC  20234 

1  AAI  Corporation 

ATTN:  H.  W.  Schuette 
PO  Box  126 

Hunt  Valley,  MD  21030-0126 

1  ABEX  Research  Center 

ATTN:  Dr.  Michael  J.  Normandin 
65  Valley  Road 
Mahwah,  NJ  07430 

2  ADPA 

ATTN:  Donna  R.  Alexander 
Bill  King 

Two  Colonial  Place,  Suite  >00 
2101  Wilson  Boulevard 
Arlington,  VA  22201-3061 

1  ARC  Professional  Services  Group 
ATTN:  Arnold  R.  Gritzke 
5501  Backlick  Road 
Springfield,  VA  22151 

1  Adelman  Associates 
ATTN:  Herbert  S.  Weintraub 
291  North  Bernardo  Avenue 
Mountain  View,  CA  94014-5205 

2  Advanced  Marine  E'  uprises 
ATTN:  James  F.  Hess 

CPT  Frederic  S.  Hering  USN  (Ret) 
1725  Jefferson  Davis  Highway 
Suite  1300 

Arlington,  VA  22202 


No.  of 

Copies  Ory*ni«ftt,inn 

1  The  Armed  Forces  Communications  and 

Electronics  Association 
ATTN:  Kirby  Lamar,  BG(Ret) 

4400  Fair  Lakes  Court 
Fairfax,  VA  22033-3899 

2  Aero  Corporation 
ATTN:  David  S.  Eccles 

Gregg  Snyder 
P.O.  Box  92957,  M4/913 
Los  Angeles,  CA  90009 

1  AFELM,  The  Rand  Corporation 
ATTN:  Library-D 
1700-Main  Street 

Santa  Monica,  CA  90406 

2  Air  Force  Wright  Aeronautical-Labs 
ATTN:  CDJ,  CPT  Jost 

CDJ,  Joseph  Faison 

Wright-Patterson  AFB,  OH  45433-6523 

1  Alliant  Computer  Company 
ATTN:  David  Micciche 
1  Monarch  Drive 
Littleton,  MA  01460 

1  Alliant  Techsystems,  Inc. 

ATTN:  Hatem  Nasr 
Systems  and  Research  Center 
3660  Technology  Drive 
P.O.  Box  1361 
Minneapolis,  MN  55418 

1  Alliant  Techsystems,  Inc. 

ATTN:  Fred  J.  Parduhn 
7225  Northland  Drive 
Brooklyn  Park,  MN  55428 

2  Alliant  Techsystems,  Inc. 

ATTN:  Raymond  H.  Burg 

Laura  C.  Dillway 
MN38-4000 

10400  Yellow  Circle  Drive 
Minnetonka,  MN  55343 


No.  of 

Copies  Org»ni«*tinn 

1  Alliston  Gas  Turbine 
Division  of  GM 
ATTN:  Michael  Swift 
PO  Box  420,  SC  S22B 
Indianapolis,  IN  46260-0420 

1  Aluminum  Company  of  America 
ATTN:  Frank  W.  Baker 
Alcoa  Technical  Center 
Alcoa  Center,  PA  15069 

1  Analysis  and  Technology 

ATTN:  RADM  Thomas  M.  Hopkins  USN 
(Ret) 

1113  Carper  Street 
McLean,  VA  22101 

1  ANSER 

ATTN:  James  W.  McNulty 
1215  Jefferson  Davis  Highway 
Arlington,  VA  22202 

1  ARC  C-500 

ATTN:  John  H.  Bucher 
Modena  Road 
Coatesville,  PA  19320 

1  Armament  Systems,  Inc. 

ATTN:  Gerard'Zeller 
P.O.  Box  158 
211  West  Bel  Air  Avenue 
Aberdeen,  MD  21001 

1  Armored  Vehicle  Technologies 
ATTN:  Coda  M.  Edwards 
PO  Box  2057 
Warren,  MI  48090 

1  ASI  Sytems,  International 

ATTN:  Dr.  Michael  Stamatelatos 
3319  Lone  Jack  Road 
Encinitas,  CA  92024 

1  Auburn  University 

Electrical  Engineering  Department 
ATTN:  Dr.  Thomas  Shumpert 
Auburn  University,  AL  36849 
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1  A.W.  Bayer  and  Associates 

ATTN:  Albert  W.  Bayer,  President 
Marina  City  Club 
4333  Admiralty  Way 
Marina  del  Rey,  CA  90202-5469 

1  Battelle  Research  Laboratory 

ATTN:  TACTEC  Library  (J.N.  Huggins) 
505  King  Avenue 
Columbus,  Ohio  43201-2693 

1  Battelle  Research  Laboratory 

Defense  and  Space  Systems  Analysis 
ATTN:  Dr.  Richard  K.  Thatcher 
505  King  Avenue 
Columbus,  Ohio  43201-2693 

1  Battelle  Research  Laboratory 
ATTN:  Bernard  J.  Tullington 
1300  N.  17th  Street,  Suite  1520 
Arlington,  VA  22209 

3  Battelle 

Edgewood  Operations 
ATTN:  Roy  Golly 

Gene  Roecker 
Robert  Jameson 
2113  Emmorton  Park  Road 
Edgewood,  MD  21040 

1  The  BDM  Corporation 
ATTN:  Edwin  J.  Dorchak 
7915  Jones  Branch  Drive 
McLean,  VA  22102-3396 

1  The  BDM  Corporation 
ATTN:  Fred  J.  Michel 
1300  N.  17th  Street 
Arlington,  VA  22209 

1  Bell  Helicopter,  Textron 
ATTNj  Jack  R.  Johnson 
PO  Box  482 
Fort  Worth,  TX  76101 


3  BMY,  Division  of  Harsco 

ATTN:  William  J.  Wagner,  Jr. 

Ronald  W.  Jenkins 
EdMagalski 
POBox  1512 
York,  PA  17404 

1  Board  on  Army  Science  and  Technology 
National  Research  Council 

Room  MH  280 

2101  Constitution  Avenue,  NW 
Washington,  DC  20418 

2  Boeing  Aerospace 

ATTN:  Dr.  Robert  Chiavetta 
Dr.  John  Kurns 
Mail  Stop  8IC17 
P"0.  Box  3999 
Seattle,  WA  98124-2499 

2  Boeing  Corporation 

ATTN:  MS  33-04,  Robert  Bristow 

MS  48-88,  Wayne  Hammond 
PO  Box  3707 
Seattle,  WA  98124-2207 

1  Boeing  Vertol  Company 
A  Division  of  Boeing  Co. 

ATTN:  MS  P30-27,  John  E.  Lyons 
PO  Box  16858 
Philadelphia,  PA  19142 

1  Booz-Allen  and  Hamilton,  Inc. 

ATTN:  Dr.  Richard  B.  Benjamin 
Suite  131,  4141  Colonel  Glenn  Hwy. 
Dayton,  OH  45431 

1  Booz-Allen  and  Hamilton,  Inc. 

ATTN:  Lee  F.  Mallett 

1300  N.  17th  Street,  Suite  1610 
Rosslyn,  VA  22209 

2  Booz-Allen  and  Hamilton,  Inc. 

ATTN:  John  M.  Vice 
WRDC/FIVS/SURVIAC 

Bldg  45,  Area  B  -  ----- 

Wright-featterson  AFB,  OH 
45433-6553 
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1  John  Brown  Associates 
ATTN:  Dr.  John  A.  Brown 
PO  Box  145 

Berkeley  Heights,  NJ  07922-0145 

1  Chamberlain 

ATTN:  Mark  A.  Sackett 
PO  Box  2545 
Waterloo,  IA  50704 

1  Commander 

Combined  Arms  Combat  Development 
ATTN:  ATZL-CAP  (LTC  Morrison) 
Dir,  Surv  Task  Force 
Ft.  Leavenworth,  KS  66027-5300 

1  Commander 

Combined  Arms  Combat  Development. 
ATTN:  ATZL-HFM:(Dwain  Skelton) 
Ft.  Leavenworth,  KS  66027-5300 

1  Computer  Sciences  Corporation 

200  Sparkman  Drive 
Huntsville,  AL  35805 

3  Computervision  Corporation 
ATTN:  A.  Bhide 

V.  Geisberg 
R.  Hillyard 

201  Burlington  Road 
Bedford,  MA  01730 

1  Cray  Research,  Inc. 

ATTN:  William  W.  Kritlow 
P.O.  Box  151 

Huntington  Beach,  CA  92648-0151 

1  CRS  Sirrine,  Inc. 

ATTN:  Dr.  James  C.  Smith 
PO  Box  22427 
1177  West  Loop  South 
Houston,  TX  77227 

1  CSC 

ATTN:  Abner  W.  Lee 
200Sparkman  Drive 
Huntsville,  AL  35805 
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2  Cypress  International 

ATTN:  August  J.  Caponecchi 
James  Logan 
1201  E.  Abingdon  Drive 
Alexandria,  VA  22314 

1  DATA  Networks,  Inc. 

ATTN:  William  E.  Regan,  Jr. 
President 

288  Greenspring  Station 
Brooklandville,  MD  21022 

J  DNA 

ATTN:  LCDR  Charles  Nofziger 
6801  Telegraph  Road 
Alexandria,  VA  22310 

1  Datatec,  Inc. 

ATTN:  Donald  E.  Gudney 
President 
326  Green  Acres 
Fort  Walton,  EL  32548 

1  David  Taylor  Research  Center 
ATTN:  Dr.  Fred  J.  Fisch 
2203  Eastlake  Road 
Timonium,  MD  21093-5000 

1  David  Taylor  Research  Center 
ATTN:  Robert  E.  Fuss 
UERD,  Code  177 
Portsmoutn,  VA  23709-5000 

1  David  Taylor  Research  Center 
ATTN:  Seymour  N.  Goldstein 
Code  1210 

Bethesda,  MD  20084-5000 

1  David  Taylor  Research'Center 
ATTN:  Ib  S.  Hansen 
Code  174 

Bethesda,  MD  20084-5000 

1  David  Taylor  Research  Center 
ATTN:  Harry  Price  Cray 
Code  1740.4 

Bethesda,  MD  20084-5000 
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1  David  Taylor  Research  Center 
ATTN:  Jackson  T.  Hawkins 
Code  1740.2 

Bethesda,  MD  20084-5000 

1  David  Taylor  Research  Center 
ATTN:  Steven  L.  Cohen 
Code  1230 

Bethesda,  MD  20084-5000 

1  David  Taylor  Research  Center 
ATTN:  Dennis  Clark 
Code  0111 

Bethesda,  MD  20084-5000 

1  David  Taylor  Research  Center 
ATTN:  John  R.  Krezel 
UERD,  Code  177.2 
Portsmouth,  VA  23709-5000 

1  David  Taylor  Research  Center 
ATTN:  Richard  E.  Metrey 
Code  01 

Bethesda,  MD  20084-5000 

1  David  Taylor  Research  Center 
ATTN:  Dr.  Paul  C.  St.  Hilaire 
Code  1210 

Bethesda,  MD  20084-5000 

1  David  Taylor  Research  Center 
ATTN:  Arthur  Marchand 
Code  2843 

Annapolis,  MD  21042 

1  David  Taylor  Research  Center 
ATTN:  Michael  Riley 
UERD,  Code  177 
Portsmouth,  VA  23709r5000 

1  David  Taylor  Research  Center 
ATTN;  J.  William  Sykes 
Code  175 

Bethesda,  MD  20084-5000 

1  David  Taylor  Research  Center 
ATTN:  Herbert  Wolk 
Code  1740.1 

Bethesda,  MD  20084-5000 
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1  University  of  Dayton 

Graduate  Engineering  and  Research 
Kettering  Lab  262 
ATTN:  Dr.  Gary  Thiele,  Director 
Dayton,  OH  45469 

1  Defense  Nuclear  Agency 
Structural  Dynamics  Section 
ATTN:  Tom  Tsai 
Washington,  DC  20305 

1  Del co  Systems  Operation 
ATTN:  John  Steen 
6767  Hollister  Avenue,  #P202 
Goleta,  CA  93117 

1  Denver  Research  Institute 

Target  Vulnerability  and  Survivabilil 
Laboratory 

ATTN:  Lawrence  G.  Ullyatt 
PO  Box  10127 
Denver,  CO  80210 

1  Denver  Research  Institute 
University  of  Denver 
ATTN:  Louis  E.  Smith 
University  Park 
Denver,  CO  80208 

1  Dow  Chemical,  U.S.A 

ATTN:  Dr.  P.  Richard  Stoesscr 
Contract  R&D 
1801  Building 
Midland,  MI  48674-1801 

1  Drexel  University 

ATTN:  Dr.  Pei  Chi  Chou 
College  of  Engineering 
Philadelphia,  PA  19104 

1  DuPont  Company  FPD 

ATTN:  Dr.  Oswald  R.  Bergmann 
B-1246,  1007  Market  Street 
Wilmington,  DE  19898 
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1  Dynamics  Analysis  and  Test  Associates 
ATTN:  Dr.  C.  Thomas  Saveli 
2231  Faraday  Ave 
Suite  103 

Carlsbad,  CA  92008 

1  E.  I.  Dupont  TED  FMC 

ATTN:  Richard  O.  Myers  Jr. 
Wilmington,  DE  19898 


1  USAETL/IAG 

ATTN:  Jim  Campbell 
Bldg  2592,  Room  S16 
Ft.  Belvoir,  VA  22060-5546 

1  FMC  Corporation 
ATTN:  Sidney  Kraus 
1105  Coleman  Ave,  Box  1201 
San  Jose,  CA  95108 


1  Eichelberger  Consulting  Company  3 

ATTN:  Dr.  Robert  Eichelberger 
President 

409  West  Catherine  Street 
Bel  Air,  MD  21014 

1  -  Electronic  Warfare  Associates,  Inc. 

ATTN:  William  V.  Ghiaramonte  5 

2071  Chain  Bridge  Road 
Vienna,  VA  22180 

1  Emprise,  Ltd. 

ATTN:  Bradshaw  Armendt,  Jr 
201  Crafton  Road 
Bel  Air,  MD  21014 

8  Environmental  Research  Institute  of  Michigan 
ATTN:  Mr.  K.  Augustyn 

Mr.  Ivozma  2 

Dr.  I.  La  Haie 

Mr.  R.  Horvath 

Mr.  Arnold 

Mr.  E.  Cobb 

Mr.  B.  Morey 

Mr.  M.  Bair 

PO  Box  8618  1 

Ann  Arbor,  MI  48107 

1  E-OIR  Measurements,  Inc. 

ATTN:  Russ  Moulton 

PO  Box  3348,  College  Station  1 

Fredericksburg,  VA  22402 

1  ERIM 

ATTN:  Stephen  R.  Stewart 
Exploitation  Applications  Department 
Image  Processing  Systems  Division 
PO  Box  8618 

Ann  Arbor,  MI  48107-8618 


FMC  Corporation 
ATTN:  Ronald  S.  Beck 
Martin  Lim 
Jacob  F.  Yacoub 
881  Martin  Avenue 
Santa  Clara,  CA  95052 

FMC  Corporation 
Advanced  Systems  Center  (ASC) 
ATTN:  Charles  A.  Millard 
Scott  L.  Langlie 
Herb  Theumer 
Walter  L.  Davidson 
J.E.  Alexander 
1300  South  Second  Street 
PO  Box  59043 
Minneapolis,  MN  55459 

FMC  Corporation 
Defense  Systems  Group 
ATTN:  Robert  Burt 

Dennis  R.  Nitschke 
1115  Coleman  Avenue 
San  Jose,  CA  95037 

FMC  Corporation 
Naval  Systems  Division  (NSD) 
ATTN:  MK-45,  Randall  Ellis 
Minneapolis,  MN  5542! 

FMC  Corporation 
Northern  Ordnance  Division 
ATTN:  M3- 11,  Barry  Brown 
4800  East  River  Road 
Minneapolis,  MN  55421 
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6  FMC  Corporation 

Ordnance  Engineering  Division 
ATTN:  H.  Croft 

M.  Hatcher 
L.  House 
J.  Jackson 
E.  Maddox 
R.  Musante 

1105  Coleman  Ave,  Box  1201 
San  Jose,  CA  95108 

1  GE  Aircraft  Engines 

ATTN:  Dr.  Roger  B.  Dunn 
One  Neumann  Way,  MD  J185 
Cincinnati,  OH  45215-6301 

1  General  Atomics 

ATTN:  Chester  J.  Everline, 

Staff  Engineer 
P.O.  Box  85608 
San  Diego,  GA  92138-5608 

1  General  Dynamics 

ATTN:  Dr.  Fred  Cleveland 
P.O.  Box  748 
Mail  Zone  5965 
Ft.  Worth,  TX  76101 

4  General  Dynamics 

ATTN:  MZ-43621 12,  Robert  Carter 
MZ-4362029,  Jim  Graciano 
MZ-4362055,  Gary  Jackman 
MZ-4362055,  Jay  A.  Lobb 
38500  Mound 

Sterling  Heights,  MI  48310 

3  General  Dynamics  Corporation 
ATTN:  MZ-2650,  Dave  Bergman 
MZ-2860,  John  Romanko 
MZ-2844,  Cynthia  Waters 
PO  Box  748 

Ft.  Worth,  TX  76101-0748 

1  General  Dynamics  Land  Systems 
ATTN:  Robert  Garter 
PO  Box  1804 
Warren,  MI  48090 
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1  General  Dynamics  Land  Systems 
ATTN:  Dr.  Pauius  Kersten 
POBox  1901 
Warren,  MI  48090 

1  General  Dynamics  Land  Systems 
ATTN:  William  M.  Mrdeza 
PO  Box  2045 
Warren,  MI  48090 

5  General  Dynamics  Land  Systems 
ATTN:  Richard  Auyer 
Otto  Renius 
N.  S.  Sridharan 
Dean  R.  Loftin 
Dr.  Phil  Lett 
PO  Box  2074 
Warren,  MI  48090-2074 

3  General  Motors  Corporation 
Research  Laboratories 
ATTN:  J.  Boyse 
J.  Joyce 
R.  Sarraga 
Warren,  MI  48090 

A  General  Motors  Corporation 
Military  Vehicles  Operations 
Combat  Vehicle  Center 
ATTN:  Dr.  John  A.  MacBain 
PO  Box  420  Mail  Code  01 
Indianapolis,  IN  46206-0420 

1  Gettysburg  College 
Box  405 

Gettysburg,  PA  17325 

1  Grumman  Aerospace  Corporation 
Research  and  Development  Center 
ATTN:  Dr.  Robert  T.  Brown, 

Senior  Research  Scientist 
Bethpage,  NY  11714 

1  GTRI-RAIL-MAD 

ATTN:  Mr.  Joe  Bradley 
CRB 577 

Atlanta,  GA  30332 
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1  Hughes  Associates 
ATTN:  J.  Thomas  Hughes 
2730  University  Blvd. 

Suite  902 

Wheaton,  MD  20902 

2  INEL/EG&G 
Engineer  Lab 
ATTN:  Ray  Berry 

M.  Marx  Hintze 
PO  Box  1625 
Idaho  Falls,  ID  83-151 

1  Interactive  Computer  Graphics  Center 
Rensselear  Polytechnic  Inst. 

ATTN:  M.  Wozny 
Troy,  NY  12181 

1  International  Development  Corporation 
ATTN:  Trevor  O.  Jones 
18400  Shelburne  Road 
Shaker  Heights,  OH  44118 

1  Intergraph 

National  Exploitation  Systems 
ATTN:  John  H.  Suter 
2051  Mercator  Drive 
Reston,  VA  22091-3413 

1  ISAT 

ATTN:  Roderick  Briggs 
1305  Duke  Street 
Alexandria,  VA  22314 

1  ITT  Defense 

ATTN:  Joseph  Conway 
1000  Wilson  Blvd. 

30th  Floor 

Arlington,  VA  22209 

1  Joint  Technical  Coordinating  Group 
ATTN:  Philip  Weinberg 
JTCG/AS5 
AIR-516J5 

■Washington,  DC  20361-5160 
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1  California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
ATTN:  D.  Lewis 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109 

1  Kaman  Sciences  Corporation 
ATTN:  Timothy  S.  Pendergrass 
600  Boulevard  South,  Suite  208 
Huntsville,  AL  35802 

1  Ketron,  Inc. 

ATTN:  Robert  S.  Bennett 
901  Dulaney  Valley  Rd  #220 
Baltimore,  MD  21204-2600 

1  Keweenaw  Research  Center 
Michigan  Technological 
University 

ATTN:  Bill  Reynolds 
Houghton,  MI  49931 

1  Lanxido  Armor  Products 
ATTN:  Dr.  Robert  A.  WolfTe 
Tralee  Industrial  Park 
Newark,  DE  19711 

2  Lincoln  Laboratory 
MIT 

ATTN:  Dr.  Robert  Shin 
Dr.  Chuck  Burt 
P.O.  Box  73 
Lexington,  MA  02173 

3  Lincoln  Laboratory 
MIT 

Surveillance  Systems  Group 
ATTN:  R.  Barnes 
G.  Knittel 
J.  Kong 
244  Wood  Street 
Lexington,  MA  02173-0073 

1  Lockheed-Corporation 
ATTN:  R/C. 'Smith 
PO  Box  551 
Burbank,  CA  91520 
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3  Lockheed-California  Company 
ATTN:  C.  A.  Burton 
R,  J.  Ricci 
M.  Steinberg 
Burbank,  CA  91520 

2  Lockheed-Georgia  Company 
ATTN:  Ottis  F.  Teuton 
J.  Tulkoff 

Dept.  72-91,  Zone  419 
Marietta,  GA  30063 

1  Lockheed  Palo  Alto  Research  Lab 
ATTN:  John  A.  DeRuntz,  JR 
0/93,  B/25I 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

1  Logistics  Management  Institute 
ATTN:  Edward  D.  Simms  Jr. 

6400  Goldsboro  Road 
Bethesda,  MD  20817-5886 

1  Los  Alamos  Technical  Associates,  Inc. 
ATTN:  JohnS.  Daly 

6501  Americas  Parkway,  #900 
Albuquerque,  NM  87110 

2  Los  Alamos  Technical  Associates,  Inc. 
ATTN:  James  C.  Jacobs 

Donald  M.  Lund 
8550  Arlington  Boulevard 
Suite  301 

Fairfax,  VA  22031 

1  Los  Alamos  Technical  Associates,  Inc. 
ATTN:  Thomas  Giacofci 
3020  Hamaker  Court 
Fairfax,  VA  22031 

1  LTV  Aerospace  and  Defense  Company 
ATTN:  Daniel  M.  Reedy 
PO  Box  655907 
Dallas,  TX  75265-5907 
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3  Martin  Marietta  Aerospace 
ATTN:  MP-113,  Dan  Dorfman 

MP-433,  Richard  S.  Dowd 
MP-243,  Thomas  C.  D’Isepo 
PO  Box  555837 
Orlando,  FL  32855-5837 

3  Mathematical  Applications  Group,  Inc. 
ATTN:  M.  Cohen 

R.  Goldstein 
H.  Steinberg 
3  Westchester  Plaza 
Elmsford,  NY  10523 

1  Maxwell  Laboratories,  Inc. 

ATTN:  Dr.  Michael  Holland 

8888  Balboa  Avenue 

San  Diego,  CA  92123-1506 

1  McDonnell  Douglas  Astronamic 
ATTN:  Nikolai  A.  Louie 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 

1  McDonnell  Douglas,  Inc. 

ATTN:  David  Hamilton 

PO  Box  516 

St.  Louis,  MO  63166 

1  McDonnell  Douglas,  Inc. 

ATTN:  Alan  R.  Parker 
3855  Lakewood  Blvd.,  MC  35-18 
Long  Beach,  CA  90846 

1  Micro  Electronics  of  North  Carolina 
ATTN:  Gershon  Kedem 
PO  Box  12889 

Research  Triangle  Park,  NC  07709 
1  MIT 

ATTN:  Dr.  S.  Benton 
RE15-416 

Cambridge,  MA  02139 
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5  The  MITRE  Corporation 

ATTN:  Edward  C.  Brady,  Vice  President 
Dr.  Robert  Henderson 
Dr.  Nicklas  Gramenopoulos 
Dr.  Narayana  Srinivasan 
Norman  W.  Huddy 
7525  Colshire  Drive 
McLean,  VA  22102-3184 

2  NFK  Engineering,  Inc. 

ATTN:  Dr.  Michael  P.  Pakstys 
Justin  W.  Held 
4200  Wilson  Blvd. 

Arlington,  VA  22203-1800 

1  NFK  Engineering,  Inc. 

ATTN:  John  J.  Turner 
1125  Trotting  Horse  Lane 
Great  Falls,  VA  22066 

1  NAS  A- Ames  Research  Center 

ATTN:  Dr.  Alex  Woo 
Mail  Stop  227-2 
Moffett  Field,  CA  94035-1000 

1  NASA-Ames  Research  Center 
ATTN:  Leroy  Presley 
Mail  Stop  227-4 
Moffett  Field,  CA  94035-1000 

1  NAVIR  DEVCON 

ATTN:  Frank  Wenograd 
Code  6043 

Walminstor,  PA  18974 

1  North  Aircraft 

ATTN:  Dr.  Athanosis  Varvatsis 
Mail  Zone  3622/84 
1  Northrop  Ave 
Hawthorne,  CA  90250 

1  Northrop  Corporation 

Research  and  Technology  Center 
ATTN:  James  R.  Reis 
One-Research  Park 
Palos  Verdes  Peninsula,  CA  90274 


1  Norton  Company 

ATTN:  Ronald  K.  Bart 
1  New  Bond  Street 
Worcester,  MA  01606-2698 

1  The  Oceanus  Company 

ATTN:  RADM  Robert  H.  Gormlev, 
(Ret) 

PO  Box  7069 
Menlo  Park.  CA  94026 

1  Oklahoma  State  University 

College  of  Engineering,  Architecture 
and  Technology 

ATTN:  Thomas  M.  Browder,  Jr. 

PO  Box  1925 

Eglin  AFB,  FL  32542 

1  Pacific  Scientific/Htl  Division 
ATTN:  Robert  F.  Aldrich 
1800  Highland  Avenue 
Duarte,  CA  91010 

1  Perceptronics,  Inc. 

ATTN:  Dean  R.  Loftin 
21111  Erwin  Street 
Woodland  Hills,  CA  91367 

1  Princeton  University 

Mathematics  Department 
Fine  Hall 
Washington  Road 
ATTN:  John  Tukey 
Princeton,  NJ  08544-1000 

1  PRI,  Inc. 

ATTN:  W.  Bushell 
Building  E4435,  Second  Floor 
Edgewood  Area-APG,  MD  21010 

1  RGB  Associates,  Inc. 

ATTN:  R.  Barakat 
Box  B 

Wayland,  MA  01778 


39 


No.  of 

Copies  Organisation 

1  Rockwell  International  Corporation 
ATTN:  Dr.  H.  Bran  Tran 
P.O.  Box  92098 
Department  113/GB01 
Los  Angeles,  CA  90009 

x  Rockwell  international  Corporation: 
ATTN:  Keith  R.  Rathjen, 

Vice  President 

3370  Miraloma  Avenue  (031-HA01) 
Anaheim,  GA  92803-3105 

1  Rome  Air  Development  Center 

ATTN:  RADC/IRRE,  Peter  J.  Costianes 
Griffis  Air  Force  Base,  NY  13441-5700 


1  RomenAir  Development  Center 

RADC/OCTM 
ATTN:  Edward  Starczewski 
Building  106 

Griffis  Air  Force  Base,  NY  13441-5700 
1  S-Cubed 

ATTN:  Michael  S.  Lancaster 
1800  Diagonal  Road,  Suite  420 
Alexandria,  VA  22314 

1  Sachs/Freeman  Associates,  Inc. 

ATTN:  Donald  W.  Lynch 

Senior  Research  Physicist 
205,  Yoakum  Parkway,  #511 
Alexandria,  VA  22304 

1  SAIC 

ATTN:  Dr.  Alan  J.  Toepfer 
2109  Air  Park  Drive,  SE 
Albuquerque,  NM  87106 

1  SAIC 

ATTN:  John  H.  McNeilly, 

.  Senior  Scientist 
171Q  Goodridge  Drive 
McLean,  VA  22102 
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2  SAIC 

ATTN:  Terry  Keller 
Robert  Turner 
Suite  200 

1010  Woodman  Drive 
Dayton,  OH  45432 

1  SAIC 

ATTN:  David  R.  Garfinkle 
Malibu  Canyon  Business  Park 
26679  W.  Agoura  Road,  Suite  200 
Calabasas,  CA  91302 

2  George-Sharp  Company 
ATTN:  Dennis  M.  McCarley 

Roger  O.  Mau 
2121  Crystal  Drive 
Suite  714 

Arlington,  VA  22202 

1  Sidwell-Ross  and  Associates,  Inc. 
ATTN:  LTG  Marion  C.  Ross, 
(USA  Ret) 

Executive  Vice  President 
PO  Box  88531 
Atlanta,  GA  30338 

1  Sigma  Research  Inc. 

ATTN:  Dr,  Richard  Boss! 
4014Hampton  Way 
Kent,  WA  98032 

1  Simula,  Inc. 

ATTN:  Joseph  W.  Coltman 
10016  South  51st  Street 
Pheonix,  AZ  85044 

1  SimTech 

ATTN:  Dr.  Annie  V.  Saylor 
3307, Bob  Wallace  Ave.,  Suite  4 
Huntsville,  AL  35807 

1  Alan  Smolen  and  Associates,  ,  Inc. 
ATTN:  Alan  Smolen,  President. 
One  Cynthia  Court 
Palm  Coast,  EL  32027-8172 
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3  Southwest  Research  Institute 
ATTN:  Martin  Goland 
Alex  B.  Wenzel 
Patrick  H.  Zabel 
P.O.  Drawer  28255 
San  Antonio,  TX  78238-0255 

3  Sparta,  Inc. 

ATTN:  David  M.  McKinley 
Robert  E.  O’Connor 
Karen  M.  Rooney 
4901  Corporate  Drive 
Huntsville,  AL  35805-6201 

1  SRI  International 

ATTN:  Donald  R.  Curran 
333  Ravenswood  Ave. 

Menlo  Park,  CA  94025 

1  Star  Laboratory,  Stanford  University 
ATTN:  Dr.  Joseph  W.  Goodman 
Electrical  Engineering  Department 
233  Durand  Building 
Stanford,  CA  94305-4055 

1  University  of  Michigan 

ATTN:  Dr.  John  F.  Vesecky 
2212  Space  Research  Blvd. 

Ann  Arbor,  MI  48109-2143 

1  Princeton  University 
ATTN:  Dr.  Curt  Callen 
Physics  Department 
PO  Box  708 
Princeton,  NJ  08544 

1  University  of  California,  San  Diego 
ATTN:  Dr.  Gordon  J.  MacDonald 
Institute  on  Global  Conflict 
and  Cooperation  (0518) 

9500  Gilman  Drive 
La  Jolla,  CA  92093-0518 


3  Structural  Dynamics  Research 
Corporation'  (SDRC) 

ATTN:  R.  Ard 

W.  McClelland 
J.  Osborn 

2000  Eastman  Drive 
Milford,  OH  45150 

1  Syracuse  Research  Group 
ATTN:  Dr.  Chung-Chi  Cha 
Merrill  Lane 
Syracuse,  NY  13210 

1  System  Planning  Corporation 
ATTN:  Ann-Hafer 
1500  Wilson  Blvd 
Arlington,  VA  22209 

1  S-Cubed 

ATTN:  Robert  T.  Sedgwick 

PO  Box  1620 

La  Jolla,  CA  92038-1620 

2  TASC 

ATTN:  Charles  E.  Clucus 
Darrell  James 
970  Mar-Wait  Drive 
Ft.  Walton  Beach,  FL  32548 

1  TASC 

ATTN:  Harry  I.  Nimon,  Jr 
1700  N.  Moore  Street,  Suite  1220 
Arlington,  VA  22209 

1  TASC 

ATTN:  COL  James  Logan  (Ret) 
1101  WilsonBlvd. 

Suite  1500 

Arlington,  VA  22209 

1  COLSA,  Inc. 

ATTN:  Mr.  Willy  Albanes 
P.O.  Box  1068 
Huntsville,  AL  35807-3301 

1  Techmatics;  Inc. 

ATTN:  Ronald  R.  Rickwald 
2231  Crystal:  Drive 
Arlington,  VA  22202-3742 
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1  Technical  Solutions,  Inc 
ATTN:  John  R.  Robbins 
P.O.Box  1148 
Mesillia  Park,  NM  88047 

1  Teledyne  Brown  Engineering 
ATTN:  John  W.  Wolfsberger,  Jr. 
Cummings  Research  Park 
300  Sparkman  Drive,  NW 
PO  Box  070007 
Huntsville,  AL  35807-7007 

1  Tradeways,  Ltd. 

ATTN:  Joseph  G.  Gorski, 

President 

307F  Mrple  Avenue  West 
Vienna,  VA  22180 

1  Ultramet 

ATTN:  Dr.  Jacob  J.  Stiglich 
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